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ABSTRACT 


Laser-induced  structural  damage  to  dielectric  thin- 
film  coatings  was  investigated  using  a  TEMQQ  Q-switv.hed  ruby  laser. 
This  study  included  mono-,  bi-,  and  multi-layered  coatings  of  the 
materials  TiC>2,  SiC>2,  ZrC>2,  MgF and  ZnS  on  substrates  of  glass, 
fused  silica,  rocksalt  and  spinel.  The  samples  included  coatings  for 
antireflection  and  reflection  at  the  ruby  wavelength  and  were  produced 
either  at  USC  or  by  a  commercial  vendor.  Damage  threshold  energy 
densities  were  measured  with  attention  paid  to  dependences  on  laser 
beam  spot-size,  film  material  and  thickness,  substrate  condition  and 
incipient  scattering  of  the  film.  Two  major  facts  about  laser  damage 
to  thin  films  were  discovered.  One  fact  is  that  laser-induced  scatter, 
indicative  of  film  break-up,  is  observed  prior  to  the  threshold  of  spark 
formation  at  the  film.  The  other  fact  is  that  the  damage  threshold  of 
the  thin-film  increases  as  the  laser  beam  spot-size  decreases.  A 
theoretical  model  was  developed  which  correlated  the  nature  and  distri¬ 
bution  of  coating  defects  to  this  spot-size  dependence. 
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I.  INTRODUCTION 


Progress  in  the  development  and  application  >f  high 
power  lasers  is  presently  limited  by  catastrophic  material  failure 
induced  by  laser  radiation.  Since  1969  an  annual  symposium, 
"Damage  in  Laser  Materials,  "  has  been  held  at  Lie  National  Bureau 
of  Standards  in  Boulder,  Colorado  and  much  effort  has  been  directed 
towards  understanding  laser -induced  damage  processes  in  transparent 
dielectrics.  Experiments  typically  use  an  intense  laser  beam,  pro¬ 
duced  by  a  high  power  solid  state  laser  with  monitored  beam  char¬ 
acteristics,  to  strike  a  small  area  of  the  tested  sample.  Measure¬ 
ments  of  damage  threshold,  in  terms  of  cither  power  density  or 
energy  density  for  different  optical  materials,  are  obtained  by  using 
various  diagnostic  techniques  for  the  detection.  The  most  frequently 
tested  specimens  are  ruby,  glass,  nonlinear  optical  crystals,  dielec¬ 
tric  thin  film  coatings  and  ir  windows.  In  solid  state  laser  hosts, 
there  has  been  substantial  progress  in  understanding  the  physics  of 
damage  and  in  most  cases,  it  can  be  reduced  to  a  problem  of 
quality  control.  However,  for  dielectric  thin-film  coatings,  the 
understanding  </f  the  damage  mechanisms  is  tentative  and  speculative. 
As  a  result,  dielectric  thin-film  coatings  are  presently  one  of  the 
weaker  components  of  high-power  laser  systems. 

The  present  work  has  been  concerned  with  the  problem 
of  developing  reliable  diagnostic  tools  to  study  laser-induced  damage 
phenomena  in  dielectric  thin  film  coatings,  and  to  understand  these 
damage  processes  in  terms  of  some  important  physical  parameters. 
The  primary  purpose  of  this  exploratory  research  was  not  so  much 
to  find  new  coating  materials  with  high  damage  thresholds,  but  to 
obtain  accurate  measurements  of  the  thresholds  of  the  best  materials 
presently  used  in  various  coating  configurations  and  with  various 
laser  irradiation  conditions.  By  examining  these  thresholds  for 
correlations  with  the  film  properties  or  laser  parameters,  it  was 


1. 


2. 


hoped  that  the  causes  of  the  damage  might  be  revealed. 

It  was  noted  in  the  laser  literature  that  a  great  amount 
of  the  past  experimental  data  were  generated  with  laser  sources  for 
which  there  was  inadequate  control  or  knowledge  of  the  laser  output, 
especially  the  transverse  spatial  profile  of  the  intensity.  At  the  early 
stage  of  this  program,  a  careful  generation  and  characterization  of 
high  power  pulses  from  a  ruby  laser  operating  in  a  single  transverse 
and  longitudinal  TEM^^  mode  was  established.  This  laser  was  used 
to  obtain  accurate  measurements,  under  well-defined  coriitions,  of 
damage  to  thin  film  coatings.  In  addition,  focusing  of  the  single¬ 
mode  pulses  also  required  special  analysis.  Equipped  with  this  well- 
controlled  high-power  tool,  we  have  directed  our  attention  to  the 
definition  of  the  threshold  of  damage  for  dielectric  thin  film  coatings 
and  several  detection  methods  were  evaluated.  The  damage  thresholds 
of  a  selection  of  vacuum  deposited  coatings  were  measured  as  a 
function  of  several  parameters  of  the  laser  beam  and  the  coating 
materials.  Many  important  and  significant  results  were  derived 
from  this  exploratory  study. 

In  Section  II,  the  experimental  setup  as  well  as  the 
laser  arrangement  with  its  characteristics  will  be  described.  The 
truncated  Gaussian  optics  related  to  the  focusing  beam  of  the 
single-mode  laser  will  be  also  discussed.  In  addition,  the  prepara¬ 
tion  of  test  specimens  which  included  mono-,  bi-,  and  multi¬ 
layered  vacuum-deposited  coatings  of  the  materials  MgF^,  SiO^, 

ZrO_,  TiOv  ZnS,  ThF  and  CeO  on  substrates  of  glass,  fused 
silica,  rock  salt,  and  spinel  will  be  presented. 

One  of  the  important  aspects  of  this  program  has  been 
the  determination  of  a  sensitive  and  reliable  criterion  for  the  onset 
of  damage.  We  have  used  four  types  of  detection  methods  for 
evaluation.  They  are:  1)  photoelectric  detection  of  the  spark; 

2)  visual  observation  of  laser-induced  increases  in  the  weak-signal 
film  scatter;  3)  microscopic  observation  of  film  breakup;  and 
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I.  INTRODUCTION 

Progress  in  the  development  and  application  of  high 
power  lasers  is  presently  limited  by  catastrophic  material  failure 
induceo  by  laser  radiation.  Since  1 9t>9  an  annual  symposium, 
"Damage  in  Laser  Materials,  "  has  been  held  at  the  National  Bureau 
of  Standards  in  Boulder,  Colorado  and  much  effort  has  been  directed 
towards  understanding  la  ser -induced  damage  processes  in  transparent 
dielectrics.  Experiments  typically  use  an  intense  laser  beam,  pro¬ 
duced  by  a  high  power  solid  state  laser  with  monitored  beam  char¬ 
acteristics,  to  strike  a  small  area  of  the  tested  sample.  Measure¬ 
ments  of  damage  threshold,  in  terms  of  either  power  density  or 
energy  density  for  different  optical  materials,  are  obtained  by  using 
various  diagnostic  techniques  for  the  detection.  The  most  frequently 
tested  specimens  are  ruby,  glass,  nonlinear  optical  crystals,  dielec¬ 
tric  thin  film  coatings  and  ir  windows.  In  solid  state  laser  hosts, 
there  has  been  substantial  progress  in  understanding  the  physics  of 
damage  and  in  most  cases,  it  can  be  reduced  to  a  problem  of 
quality  control.  However,  for  dielectric  thin-film  coatings,  the 
understanding  of  the  damage  mechanisms  is  tentative  and  speculative. 
As  a  result,  dielectric  thin-film  coatings  are  presently  one  of  the 
weaker  components  of  high-pewer  laser  systems. 

The  present  work  has  been  concerned  with  the  problem 
of  developing  reliable  diagnostic  tools  to  study  laser-induced  damage 
phenomena  in  dielectric  thin  film  coatings,  and  to  understand  these 
damage  processes  in  terms  of  some  important  physical  parameters. 
I’he  primary  purpose  of  this  exploratory  research  was  not  so  much 
to  find  new  coating  materials  with  high  damage  thresholds,  but  to 
obtain  accurate  measurements  of  the  thresholds  of  the  best  materials 
presently  used  in  various  coating  configurations  and  with  various 
laser  irradiation  conditions.  By  examining  these  thresholds  for 
correlations  with  the  film  properties  or  laser  parameters,  it  was 
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4)  time-dependent  observation  of  laser -induced  damage  phenomena. 

In  Section  III,  a  detailed  description  of  these  methods  and  their  eval¬ 
uation  will  be  given.  Light  scattering  induced  by  a  low  intensity 
laser  (e.  g.  ,  a  He-Ne  laser)  was  the  most  sensitive  sign  of  film 
breakup  for  coatings  of  low  initial  scatter.  A  correlation  between 
the  damage  thresholds  and  diffuse  weak-signal  scattering  of  multilayer 
reflectors  was  also  established.  Time-resolved  measurements  of 
scatter  have  also  shown  that  the  thin-film  damage  occurred  within 
the  laser  pulsewidth  of  20  nsec. 

In  the  course  of  'this  study,  the  importance  of  the 
electric  field  distributions  of  the  tested  specimens  was  demonstrated. 

A  computer  program  was  set  up  to  calculate  the  electric  field  dis¬ 
tribution  for  all  the  thin  film  saimples  that  were  tested  as  discussed 
in  Section  IV.  The  damage  thresholds  of  thin  films  were  strongly 
dependent  on  the  standing-wave  patterns  of  the  internal  electric  fields. 
The  entrance-face  thresholds  were  equal  to  or  greater  than  the  exit- 
face  thresholds  of  thin  films.  Furthermore,  when  the  laser  beam  is 
incident  at  an  oblique  angle,  calculation  of  field  distribution  indicates 
that  the  power-density  for  the  polarizations  parallel  and  perpendicular 
to  the  plane  of  incidence  are  substantially  different. 

In  addition  to  the  study  of  damage  to  entrance  and 
exit-face  coatings,  the  damage  threshold  energy  and  power  densities 
were  measured  as  a  function  of:  (1)  laser  beam  spot -size,  (2)  pulse- 
width  of  laser  beam  (3)  transverse  mode  structure  of  laser  beam, 

(4)  coating  material,  (5)  single  layer  thickness  of  coating,  (6)  sub¬ 
strate  material  and,  (7)  multilayer  configuration.  One  of  the  most 
important  results  in  this  study  is  that  the  damage  threshold  of  the 
thin  film  increases  as  the  spot-size  of  the  laser  beam  decreases. 

The  spot-size  is  the  1/e^  radius  of  the  intensity  profile  of  the  TEM^ 
ruby  laser  beam.  A  simple  model  has  been  developed  in  correlating 
the  nature  and  distribution  of  coating  defects  to  this  spot-size  depen¬ 
dence,  e.  g. ,  the  probability  of  the  laser  beam  striking  a  defect  site 
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will  be  greater  for  larger  spot-sizes,  while  damage  in  materials  can 
be  distinguished  as  defect  damage  and  intrinsic  damage.  The  detailed 
account  and  experimental  confirmation  for  this  model  can  be  found  in 
Section  V.  Experimental  results  of  other  laser  parameters  discussed 
in  Section  VI  are  tentative. 

The  experimental  results  of  damage  thresholds  on 
varying  several  coating  parameters  are  given  in  Section  VII.  In 
testing  single  layer  thin  films,  the  damage  thresholds  for  films  with 
low  refractive  indices  were  greater  than  those  for  films  with  high 
indices  at  the  ruby  laser  wavelength.  The  role  of  absorption  as  a 
damage  mechanism  was  analyzed.  It  was  suggested  that  linear 
absorption  could  raise  the  film  temperature  either  to  the  melting 
point  or  to  a  value  high  enough  to  produce  harmful  thermostrictive 
forces.  Section  VIII  summarizes  some  important  findings  of  this 
program  arid  some  suggestions  for  further  work  in  this  research 
area  are  mad;,. 
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II.  EXPERIMENTAL  PROCEDURES 

1.  Experimental  Setup 

In  order  to  accurately  interpret  the  results  of  experi¬ 
ments  involving  laser-induced  damage,  the  temporal  and  spatial- 
intensity  profiles  of  the  laser  pulse  must  be  smooth  and  preferably 
Gaussian.  Accordingly,  a  high  power,  passively  Q-switched  ruby 
laser  oscillator  and  amplifier  system  was  constructed  which  operated 
in  a  single  longitudinal  and  transverse  TEMQ0  mode.  The  time 
contour  of  the  pulse  was  nearly  Gaussian  with  pulsewidths  (FWHM) 
from  9  to  35  nsec,  and  peak  powers  up  to  10  MW  when  using  the 
amplifier.  The  output  characteristics  of  this  laser  system  are  listed 
in  Table  1.  The  values  given  for  the  oscillator  were  determined  at 
100  cm  (1.1  Rayleigh  distances)  from  the  oscillator  aperture.  For 
such  Gaussian  pulses,  the  peak  power  Pq  is  related  to  the  total  energy 
E  by  Pq  =  (0. 941E)/T,  where  T  is  the  full  width  of  the  pulse  at  half 
maximum.  Absolute  energy  calibration  was  performed  using  a 
ballistic  thermopile  (TRG  100),  and  agreed  to  within  5%  of  the  calibra¬ 
tion  of  another  thermopile  used  by  C.R.  Guiliano  at  Hughes  Research 
Laboratories.  In  the  far  field  of  the  laser  oscillator,  the  spatial 
intensity  profile  of  the  beam  was  close  to  a  Gaussian  distribution, 
as  determined  by  a  pinhole  scan.  At  the  distance  of  100  cm 
(1.1  Rayleigh  distance)  from  the  oscillator,  the  intensity  distribution 
at  the  temporal  peak  of  the  pulse  is  shown  in  Figure  1.  The 
amplif  ed  profile  also  had  a  similar  near-Gaussian  shape  as  shown 
in  Figure  2. 

For  measuring  damage  thresholds  of  dielectric  thin 
film  coatings,  the  laser  was  arranged  as  shown  in  Figure  3.  Most 
of  the  tests  were  performed  with  the  output  of  the  oscillator  alone, 
since  its  spatial  profile  in  the  far  field  was  near  ideal  for  analytical 
computations.  For  time  resolved  measurements,  another  laser  having 
similar  output  characteristics  but  with  only  the  oscillator  was  setup 
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NORMALIZED  SPATIAL  INTENSITY  DISTRIBUTION 
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Fiture  2.  Normalized,  transverse  spatial-intensity  distribution  of  the  ruby  amplifier.  The 
amplifier  input  face  was  at  2.2  R  (200  cm)  from  the  oscillator. 
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Figure  3.  Schematic  of  the  Q-switched  ruby  laser  system  and  spark  detection 
electronics  for  damage  threshold  measurements. 
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as  shown  in  Figure  4.  However,  unless  specified  otherwise,  our 
discussion  will  be  mainly  on  the  first  experimental  setup.  The 
detailed  operation  of  this  apparatus  was  discussed  at  the  Third 
ASTM-NBS  Damage  Symposium  [  l]. 
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Table  1.  Single-mode  output  of  the  mby  laser. 


Parameter 

Oscillator 

Amplifier 

Energy 

5-7 

mJ  (  -  5%) 

up  to  120  mJ 

Puisewidth 

** 

12  nsec 

~  12  nsec 

Peak  powe  r 

0.4  - 

0.6  MW 

<  10  MW 

Divergence  beam 

1.  1  mrad  (full  angle) 

Beam  spot- size 

0.62 

-  0.  02  mm 

Peak  energy  density 

0.  8 

-  1.  2  J/cm2 

< 

100  J/cm2 

Peak  power  density 

65  - 

100  MW/cm2 

< 

10  GW/cm2 

2.  Calibration  and  'f  uncatcd  Gaussian  Optics 
The  spatial  intensity  distribution  of  the  laser  beam 
cannot  be  exactly  Gaussian  since  the  single  mode  is  coupled  out  of 
the  laser  oscillator  through  a  circular  aperture,  resulting  in  the 
truncation  of  the  aperture  intensity  illumination.  This  truncation 
produces  rings  in  the  far-field  pattern  and  structure  within  the 
central  disk  in  the  near -field.  Figure  5  shows  the  diffraction 
pattern  of  the  laser  oscillator  in  the  far-field  with  three  different 
exposures.  The  central  disk  is  nearly  Gaussian  and  in  the  over- 
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Photographs  of  the  diffraction  pattern  of  the  giant-pulse  ruby 
laser  oscillator  (N  =  0.4)  at  3.0  Rayleigh  distances  (far  field) 
from  the  oscillator  aperture.  The  three  photographs  were 
exposed  tinder  three  different  beam  attenuations. 
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exposed  pictures,  the  presence  of  the  rings  around  the  central  disk 

due  to  truncation  is  strikingly  evident.  Now,  the  importance  of 

recognizing  the  existence  of  truncation  effects  lies  in  the  fact  that, 

if  attention  is  not  paid  to  these  effects,  calculations  of  the  axial 

intensity  may  be  in  error  by  as  much  as  100%.  Such  discrepancies 

are  particularly  possible  when  truncated  Gaussian  beams  are  focused. 

When  a  laser  beam  is  focused  by  a  lens  the  focused 

intensities  can  be  -simply  predicted  only  if  the  spatial  profile  of  the 

beam  has  the  ideal  Gaussian  distribution  [2].  Since  an  exactly 

Gaussian  beam  remains  Gaussian  through  out  an  optical  system,  only 

the  calculation  of  the  spot- size  of  the  Gaussian  beam  with  respect  to 

2 

the  lens  is  required.  The  spot-size  is  the  1/e  radius  of  the  inten¬ 
sity  profile  of  a  beam.  The  variation  of  the  spot-size  as  the  beam 
propagates  along  the  axis  of  the  lens  is  controlled  by  the  spot-size 

w  at  the  waist  of  the  input  beam  and  the  distance  d,  of  the  input 
o  1 

waist  from  the  lens  as  well  as  the  geometric  focal  length  F  of  the 
lens.  The  geometry  of  the  input  beam  to  the  lens  is  illustrated  in 
Figure  6,  and  the  value  of  the  spot-size  w^  at  a  distance  d^  from 
the  lens  is  determined  from  Eq.  (1). 
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The  parameter  b  is  defined  as  /X  where  X  is  the 

laser  wavelength.  This  Eq.  (1)  is  the  general  expression  giving 

w_  as  a  function  of  w  ,  d,  and  F,  and  can  be  used  to  determine 
2  o  1 

4  e  peak  power  and  energy  densities  anywhere  in  a  focused  beam 
provided  the  beam  has  an  ideal  Gaussian  profile.  Many  experimenters 
using  single  transverse-mode  lasers  generally  assume  that  this 
Gaussian  description  of  focused  beams  applies  to  their  laser  setup. 
Truncation  of  the  Gaussian  beam  can  severely  alter  this  calculation, 
though,  and  prevent  a  proper  analysis  of  an  experiment. 
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Truncation  of  the  Gaussian  distribution,  which  is 
infinite  in  extent,  can  never,  of  course,  be  avoided  but  is  generally 
unimportant  in  low-power  lasers.  In  high-power  Q- switched  lasers 
though,  due  to  their  resonator  design  requirements,  the  truncation 
is  usually  not  negligible.  Figure  7  shows  several  possible  intensity 
diet*  '.utions  in  t»;r  aperture  plane  of  a  laser  oscillator  having  a 
Fresnel  number  of  0.  4,  the  Fresnel  number  of  our  oscillator.  The 
profile  for  the  fundamental  stationary  (Fox-Li)  mode  is  a  truncated 
Gaussian  having  a  truncation  parameter  of  0.  84.  The  truncation 
parameter  is  the  ratio  of  the  radius  of  the  truncating  aperture  a 
to  the  spot-size  w.  This  Fox-Li  profile  is  flattened  by  saturation 
of  the  gain  [3],  which  usually  occurs  in  high-power  pulsed  lasers. 
Three  examples  are  shown  for  initial  numerical  gains  of  20,  12.2 
and  7.4,  all  calculated  assuming  a  uniform  initial  gain  distribution. 
When  the  initial  gain  distribution  is  peaked  symmetrically  about  the 
resonator  axis  due  to  the  design  of  a  particular  pumping  scheme, 
the  Fox-Li  profile  is  sharpened,  increasing  the  truncation  parameter. 
A  typical  example  of  axisymmetric  pumping  is  the  elliptical  focusing 
of  pump  light  into  a  cylindrical  laser  rod.  For  the  laser  oscillator 
used  in  our  investigations,  the  truncation  parameter  was  near  0.6. 

Since  Campbell  and  DeShazer  [4]  have  shown  that 
diffraction  effects  due  to  truncation  cannot  be  neglected  when  the 
truncation  parameter  is  less  than  2. 0,  the  usual  Gaussian  description 
of  focused  beams  does  not  apply  to  our  laser  experiment.  There¬ 
fore,  a  calculation  of  focused  truncated  Gaussian  beams  is  required 
for  the  analysis.  A  formula  has  been  published  [5j  for  the  intensity 
variation  along  the  optical  axis  for  a  special  case  of  a  truncated 
Gaussian  beam  where  the  lens  is  also  the  aperturing  element. 

U.  O.  Farrukh  [6]  has  calculated  the  diffraction  effects  of  focused 
truncated  Gaussian  beams  for  the  more  usual  situation  shown  in 
Figure  8.  The  source  laser  is  represented  by  a  collimated  ideal 
Gaussian  beam  of  spot- size  w,  apertured  by  a  diverging  lens  of 
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geometric  focal  length  and  the  radius  a.  A  lens  of  geometric 
focal  length  is  placed  a  distance  z  from  the  lens  fj ,  and  the  beam 
is  viewed  a  distance  y  away  from  the  lens.  The  axial  intensity  I  at 
the  distance  y  is  determined  from  Eq.  (2). 


Even  though  this  is  a  rather  complex  expression,  the  axial  intensities 
can  be  easily  plotted  by  computer  and  the  plots  have  predicted  the 
experimental  results. 

In  our  experiment,  a  lens  of  20.  7  cm  focal  length 
was  located  109.5  cm  from  the  laser.  The  effective  value  of  f 
was  calculated  to  be  -154.5  cm  from  the  theoretical  radial  phase 
distribution  [ 3j.  Since  the  equivalent  Gaussian  (corresponding  to 

the  measured  laser  divergence)  and  truncated  Gaussian  calculations 
of  the  axial  intensity  differ  greatly  near  the  focus,  the  axial  intensity 
was  directly  measured  and  compared  to  the  two  calculations.  The 
axial  intensity  of  such  small  spots  can  be  accurately  determined  by 
measurements  of  the  energy  required  to  damage  Polaroid  film. 

A  microscope  is  required  to  examine  the  film.  Fresh  undeveloped 
Polaroid  film  (Type  410)  was  calibrated  to  have  a  damage  threshold 
at  50  mJ/cm  ,  which  was  independent  of  the  laser  beam  spot-size. 
Figure  9  shows  the  measured  axial  intensity  distribution  after  the 
lens  and  its  comparison  to  both  the  equivalent  Gaussian  and  truncated 
Gaussian  laser  beam  calculations.  The  measured  intensity  value  near 
the  focus  is  about  twice  that  predicted  by  the  equivalent  Gaussian 
calculation.  The  difference  between  the  measured  distribution  and 
the  truncated  Gaussian  calculation  is  not  as  large  as  indicated 
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because  the  actual  truncation  parameter  is  near  0.6  and  not  the 
0.  5  as  plotted. 

Observation  of  the  intensity  profile  of  a  truncated 
Gaussian  beam  after  the  focus  reveals  axial  maxima  and  minima, 
unlike  the  smooth  distribution  before  the  focus.  This  is  explained 
by  recognizing  that  near-field  axial  extrema  are  present  before  the 
lens,  which,  when  the  lens  is  used  in  the  far-field  ( z>  2a  /X),  are 
transformed  to  locations  after  the  focus.  Figure  10  illustrated  that 
the  intensity  distribution  is  not  symmetrical  about  the  focal  plane 
for  truncation  parameters  less  than  2.  Therefore,  it  is  quite  impor¬ 
tant  to  place  the  test  samples  before  the  focus  in  order  to  know  the 
laser  power  densities  incident  to  the  samples. 

In  determining  the  damage  thresholds,  areas  of  a 
sample  were  irradiated  by  one  shot  only.  Whether  damage  occurred 
or  not,  a  new  location  was  irradiated  on  each  shot.  The  auxiliary 
He-Ne  laser  beam,  which  travelled  collinearly  with  the  pulsed  laser 
beam,  indicated  the  prospective  irradiation  sites. 

The  total  energy  and  pulsewidth  of  each  laser  pulse 
were  measured,  and  the  presence  of  a  spark  at  the  sample  surface 
was  detected  photoelectrically. 


3.  Thin  Film  Preparation 

This  study  included  mono-,  bi-,  and  multi-layered 
vacuum  deposited  coatings  of  the  materials  MgF  ,  SiO  ,  ZrO  , 

M  W 

TiO^,  ZnS,  ThF^,  and  CeO^  on  substrates  of  glass,  .’used  silica, 
rock  salt,  and  spinel.  These  samples  were  prepared  by  an  in- 
house  facility  as  well  as  several  commercial  sources.  Methods 
used  in  the  preparation  were  the  resistive-heating  evaporation, 
electron-gun  evaporation,  as  well  as  rf  sputtering.  me  thin-film 
coatings  available  for  our  experiments  are  listed  in  Tal  le  2. 
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HI.  DETECTION  METHODS 

An  important  aspect  of  this  program  has  been  the 
determination  of  a  sensitive  and  reliable  criterion  for  the  onset 
of  damage.  Many  investigators  of  laser-induced  surface  damage  to 
transparent  dielectrics  have  observed  a  luminous  surface  spark 
whenever  physical  damage  occurred.  For  example,  see  Figure  11. 
Some  have  proposed  that  the  spark  is  responsible  for  the  observed 
damage.  Although  it  has  become  a  convenient  practice  to  identify 
the  spark  threshold  as  the  damage  threshold,  other  researchers 
also  detected  damage  without  seeing  any  sparks.  Four  different 
kinds  of  detection  methods  were  studied  in  this  program:  1)  photo¬ 
electric  detection  of  the  spark,  namely,  spark  detection  system, 

2)  visual  observation  of  laser-induced  increases  in  the  weak-signal 
film  scatter,  namely,  laser-induced  scatter  (LIS),  3)  microscopic 
observation  of  film  breakup  and  damage  morphology  including  optical 
and  scanning-electron  microscopy,  4)  time  evolution  of  laser  induced 
damage,  namely,  time  resolution  optical  probe  technique.  In  the 
following,  these  methods  will  be  described  and  then  evaluated 
accordingly. 

1.  Spark  Detection  System 

The  spark  emission  was  detected  with  a  photomultiplier 
tube  (RCA  6199)  at  4340A  using  a  narrow-band  interference  filter 
(FWHM  =  40A).  Detection  of  only  the  blue  emission  decreased  the 
interference  from  the  flash-lamp  light.  Thr  photoelectric  signal 
caused  by  the  spark  emission  was  measured  with  a  Tektronix  Model 
555  oscilloscope  and  was  effectively  isolated  from  the  blue  emission 
of  the  lamp  by  means  of  a  passive  high  pass  filter.  The  electronic 
circuitry  of  the  detection  system  is  shown  in  Figure  12. 

Examples  of  oscillograms  corresponding  to  laser  - 
induced  sparks  are  illustrated  in  Figure  13.  Signals  for  the 
incident  laser  energies  causing  the  sparks  are  also  shown.  The 
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Figure  11.  Laser-induced  sparks  above  the  damage  threshold 
of  a  multilayer  reflector.  The  two  smaller  images 
shown  in  the  front  view  (a)  are  a  tiny  spark  (above) 
due  to  energy  inter -reflected  in  the  focusing  lens 
and  focused  flashlamp  illumination  (behind;  flat 
image)  reflected  to  the  camera. 


Detector:  RCA  6199  Photomultiplier 

C„B  67  pfd. ,  Rieetime  *  2. 5  nsec. 

Connecting  Cable:  RG  58/U  ,C,SI50  pfd. 

Low  Frequency  Filter;  R,sl50Kfl 

C2  *  ICO  pfd. 

Oscilloscope:  Type  IAI  Tektronix  Plug-  In 

Type  555Dual  Beam  Unit 
R2=  I  M fi  ,  C3=  15 pfd. 


Figure  1 2.  Electronic  circuit  for  photoelectronic  detection 
of  laser-induced  spark  emission. 


Figure  13.  Measurement  of  laser-induced  spark  emission 
(a)  above  and  (b)  near  the  spark  threshold. 
Also  measured  is  the  incident  laser  energy 
(upper  traces).  Dual-beam  oscillograms 
with  lOmv/cm  vertical  scale  and  time  scales 
of  200M  sec /cm  (energy)  and  50  psec/cm 
delayed  by  1  msec  (spark)  are  shown. 
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height  of  the  signal  deflection  was  a  measure  of  the  total  energy 

contained  in  the  spark  or  laser  beam.  The  sensitivity  of  the 

detection  system  is  evaluated  in  detail  in  Technical  Report  No.  5 

of  this  program.  For  a  typical  damage  spot  of  lOpim-radius  area, 

it  is  estimated  a  spark  emission  per  unit  area  and  wavelength  of 
2 

1.6  watts/cm  (am  can  be  detected.  The  temperature  corresponding 
to  this  value  is  23Q0°K,  which  is  20C°K  higher  than  the  melting 
point  of  TiO^. 


2.  Laser-Induced  Scatter  (LIS) 

Evidence  of  laser-induced  disruption  of  the  surface  of 
a  coating  was  manifested  by  increased  scattering  of  He-Ne  laser, 
travelling  collinearly  with  the  pulsed  laser  beam,  illuminated  the 
test  area.  By  visually  comparing  the  level  of  scattering  in  a 
darkened  room  before  and  after  pulsed  irradiation,  changes  in  the 
coating  structure  could  be  detected.  These  changes  were  considered 
as  damage.  For  single-and  bi-layer  coatings,  which  had  very  little 
incipient  scattering,  this  method  provided  an  especially  sensitive 
measure  of  the  damage  threshold. 

Another  approach  of  the  weak-signal  scatter  was 
measured  by  the  intensity  of  light  scattered  at  about  135°  from 
the  direction  of  the  incident  beam  using  a  severely  attenuated  output 
from  the  TEM-.  Q-switched  ruby  laser  instead  of  a  He-Ne  gas  laser. 
This  beam  had  a  power  density  of  40  MW/cm  and  did  not  produce 
any  changes  in  the  films. 

3.  Microscopy 

The  coatings  were  visually  examined  after  most  shots 
with  low-power  microscopes  (7X  or  40X).  Since  this  method  was 
much  less  sensitive  than  the  first  two  methods,  information  obtained 
was  limited  to  the  shape  and  extent  of  the  damaged  areas.  More 
detailed  information  on  the  structure  of  damage  near  threshold  was 
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provided  by  photographs  of  damaged  areas  under  high  magnification 
( 1 00X  to  800X)  with  a  Bausch  and  Lomb  Metallograph,  and  even  higher 
magnification  (up  to  20,  000)  with  a  scanning  electron  microscope  (SEM). 

4.  Time  Resolution  Technique 

The  time  evolution  of  laser-induced  damage  to  thin-films 
was  measured  by  using  an  optical  probe  technique,  as  well  as  observ¬ 
ing  the  transmitted  pulse.  This  technique  consists  of  illuminating  the 
surface  site  to  be  damaged  with  a  He-Ne  laser,  and  then  detecting 
the  reflected  beam  with  a  fast  electro-optical  system.  Any  change 
at  the  surface  of  the  sample  is  recorded  as  a  change  in  the  intensity 
of  the  reflected  probe  beam,  and  the  temporal  behavior  of  the  reflected 
probe  beam  is  directly  related  to  the  temporal  change  at  the  sample 
surface.  The  damage  is  caused  by  a  TEM^  Q-switched  ruby  pulse 
of  pulsewidth  10-20  nsec.  The  temporal  shape  of  the  ruby  pulse 
(optically  delayed)  and  the  transmitted  portion  of  the  pulse  after 
passing  through  the  sample  are  recorded  as  input  and  output  pulses 
respectively.  Each  damage  site  is  indexed  to  allow  comparison 
of  the  SEM  photograph  and  the  corresponding  time  resolved  damage 
trace.  Detailed  description  of  the  method  can  be  found  in  Quarterly 
Technical  Reports  No.  4  and  No.  7  as  well  as  reference  [7]. 

5.  Evaluation 

a.  LIS -Spark  Ratio 

It  was  found  that  for  a  single  or  bilayer  film  the  laser- 
induced  increase  in  the  film  scattering  (LIS)  was  the  most  sensitive 
method  for  detecting  damage.  (See  Table  3).  For  multilayer 
reflecting  films,  for  which  the  initial  He-Ne  scatter  level  was  quite 
high,  spark  detection  was  not  sensitive.  An  interesting  result  was 
that,  for  single  or  bilayer  films,  an  increase  in  the  He-Ne  light 
scatter  was  observed  prior  to  or  at  the  detectable  spark  threshold, 
depending  on  the  film  material.  If  the  initial  He-Ne  scatter  levels 
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had  been  low  enough  to  allow  detection  of  small  increases  in  the 
scatter  level,  the  same  result  would  have  been  expected  for  the  multi¬ 
layer  films  also. 

Table  3.  Comparison  of  spark  and  laser-induced  scattering  thresholds. 
(Laser  pulsewidths  were  nominally  12  nsec.) 


Sample 

Description 

Peak  energy 
Spark 

density 

LIS 

Spark-to-LIS 

Ratio 

J/cm^ 

J/cm^ 

U203 

ZnS,  \l  4 

27 

6 

4.  5 

0137 

TiOz,3X/4 

15 

7.  5-10.  5 

1.  5-2 

0129 

ZrOz,  X/4 

26-29 

15-19 

1.  5-2 

0133 

MgF2,  X/4 

44-53 

44-53 

1 

0135 

MgF2,  X/2 

33-40 

19-40 

1-1.  5 

0125 

Si02,  X/4 

44-51 

44-51 

1 

0127 

Si02,  X/2 

46 

46 

1 

b.  Weak-Signal  Scatter 

A  correlation  between  the  damage  thresholds  and 
diffuse  weak-signal  scattering  of  multilayer  reflectors  is  apparent 
in  Table  4.  That  the  damage  thresholds  are  generally  lower  for 
coatings  with  high  scatter  is  a  reasonable  result.  The  effective 
absorption  coefficient  of  a  coating  is  enhanced  via  internal  reflections 
of  the  scattered  energy.  If  the  contribution  of  scattering  to  the 
net  deposition  of  energy  is  great  enough,  the  resultant  damage  threshold 
would  be  lower  than  the  threshold  for  no  scattering. 


Table  4.  Comparison  of  weak-signal  scatter  and  spark  t  ireshold  for 
multilayer  reflcttors. 


Sample 

Reficctor 

Normalized  scatter 
intensities  at  94 3 A 

Peak  energy  density 
(spark) 

J/cm^ 

0104 

Ti<ysio2 

1.0 

107-127 

0102 

II 

2.  3 

98-110 

aum 

ZnS/ThF . 

4 

2.  8 

25-26 

0101 

Ti°2/Si°2 

3.  5 

83-103 

0103 

•  1 

3.  7 

121 -126 

SI  01 

It 

6.  5 

44-56 

0105 

Z'°2/Sl°2 

9.  3 

18-20 

SSI  03 

Ce02/S10. 

14.  3 

17-19 

0106 

zr02/r»io2 

16.7 

18.  5 

SS101 

CeOz/SiO? 

24.  5 

11.  5 

*0108 

Zr02/MgF2 

24.  5 

90 

*0107 

It 

27.0 

81-113 

SSI  02 

o<ys.o2 

37.  5 

14. 5-17 

SI  03 

Zr02/SI02 

41.  5 

7.  5 

*Notable  exceptions  to  the  correlation. 
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c.  Damage  Morphology 

The  visual  examination  of  irradiated  samples  using 
low-power  microscopes  (7x  and  40x)  was  not  as  sensitive  as  the  other 
methods.  However,  certain  structures  in  the  damage  morphology 
can  be  observed  by  a  high-power  optical  microscope  (lOOx  to  lOOOx) 
or  electron  microscope  (up  to  lOOKx).  In  general,  the  damaged 
areas  below  spark  threshold  were  typically  circular,  nearly  following 
the  intensity  profile  of  the  incident  beam.  Figure  14  to  19  illustrates 
that  damage  with  a  spark  shows  a  very  round  hole  of  rather  large  size, 
while  Figure  20  and  21  show  that  damage  occurs  without  always 
producing  a  spark. 


d.  Optical  Probe  Technique 

Results  presented  here  are  representative  of  entrance 
surface  damage  to  X/4  ZnS  films  on  BSC2  glass.  Damage  was  time 
resolved  with  a  resolution  of  2  nsec.  The  damage  evolution  falls 
into  two  classes:  a  fast  damage  initiation  process  starting  near  the 
peak  of  the  ruby  pulse  and  having  a  rise  time  of  ~4  nsec,  and  a 
slow  process  initiating  at  the  trailing  edge  or  after  the  ruby  pulse 
and  having  a  rise  time  10-20  nsec.  When  th-?  probe  intensity  recovers 
under  conditions  described  in  above,  the  recovery  time  is  25-50  nsec. 
Thin  film  damage  areas  on  the  order  of  4pm  in  diameter  have  been 
observed  without  spark  emission,  and  without  distortion  occurring  in 
the  transmitted  ruby  pulse.  A  typical  result  is  displayed  in  Figure  22. 
The  probe  response  indicates  a  fast  damage  initiation  process 
starting  near  the  peak  of  the  ruby  pulse  resulting  in  a  fast  decrease 
of  the  probe  intensity  to  a  minimum  value  typically  in  less  than 
4  nsec. 


Laser-induced  damage  above  the  spark  threshold 
for  a  22  layer  TiC^/SiC^  reflector  of  configurate 
G(HL)*®HL^A.  (by  Bausch  &  Lomb  Metallograph) 


. 


Laser-induced  damage  far  above  the  spark  threshold 
for  the  same  reflector  a6  in  Fig.  14. 


Severe  laser-induced  damage  on  the  same  reflector 
shown  in  Fig.  14.  The  glass  substrate  is  seen  on 
the  site,  (magnification:  llOOx) 


Laser-induced  damage  far  above  the  spark  threshold 
on  a  Ti02/Si02  multilayer  reflector  (Sample  SI  01). 
Random  penetrations  to  deeper  layers  is  apparent, 
(magnification:  120x) 
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Figure  18.  Scanning  Electron  Microscope  (SEM)  photograph  of  a 
damage  site  on  a  single  quarter-wave  thick  film  of 
Ti02  (Sample  0121)  where  a  spark  was  detected. 


Figure  19. 


SEM  photograph  of  a  damage  site  on  MgF2,  */4  film 
(Sample  0133).  The  site  dimension  is  100  microns 
from  left  to  right  (80°  sideview  of  the  site). 


Figure  21.  SEM  photograph  of  a  damage  site  on  a  single  quarter- 
wave  film  of  ZnS  (U5a)  using  US  diagnosis. 


torn* 


H 

/6  »£tc 


/Onsec 


Figure  22.  Entrance  surface  damage  to  X/4  ZnS.  (a)  Probe  trace  (upper 

trace)  and  the  damaging  ruby  pulse  (lower  trace).  The  maximum 
decrease  in  the  probe  intensity  is  35%  of  peak  intensity  (peak 
intensity  is  38  mV),  (b)  The  output  damaging  pulse  is  shown 
followed  by  the  input  delayed  reference  pulse,  (c)  SEM  photo¬ 
graph  taken  at  10°  tilt. 
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IV.  ELECTRIC  FIELDS  IN  THIN  FILM  COATINGS 


1.  Introduction 

In  the  analysis  of  laser-ind  iced  damage  processes,  the 
role  of  absorption  cannot  be  neglected.  A  more  detailed  discussion 
of  absorption  will  be  presented  in  Section  VII.  In  this  section,  we 
will  discuss  the  role  of  electric  fields  which  is  indirectly  related 
to  the  absorption  in  the  laser-induced  damage. 

The  power  P  absorbed  per  unit  volume  from  an  electro 
magnetic  wave  of  intensity  I  traveling  in  a  medium  having  a  complex 
refractive  index  N  is  P  =  -dl/dz  =  01.  0  is  the  absorption  coeffi¬ 
cient  which  is  related  to  the  imaginary  part  of  N=  n-ik  Dy  the 

—  2  —  2 

relation  0  =  4nk/X.  For  a  plane  wave  1=  E  / T)  where  E  is 
the  mean  square  of  the  electric  field  averaged  over  several  cycles 
(1  (El  ^  )  and  T|  =  (|i/  €  )  is  the  wave  impedance.  In  terms  of 
the  incident  intensity  I  (  =  |eJ  /  2 T^),  the  power  absorbed 
(watts/cm3)  at  a  distance  z  in  a  weakly  absorbing  medium  is 


P(z) 


0  n 


EM 

_+ 


I 

o 


4tt  nk 

X 


EM 


o 


2 


I 

o 


(3) 


Similar  expression  can  be  obtained  for  the  energy  absorbed  per 
unit  volume  (joules/cm  )  from  a  light  pulse  with  incident  energy 
density.  Since  the  energy  absorbed  per  unit  volume  is  proportional 
to  the  square  of  the  electric  field,  the  standing -wave  patterns  due 
to  the  reflections  at  the  film  interfaces  must  be  taken  into  account. 

2.  Calculation  of  Electric  Fields  in  Thin  Film  Coatings 

The  electric  field  intensities  were  calculated  assuming 
plane  waves  and  planar  boundaries.  In  addition,  all  films  and  sub¬ 
strates  were  assumed  to  have  negligible  absorption.  The  matrix 
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method  t8]  involving  Fresnel  coefficients  was  used  for  the  calculations. 
In  matrix  form,  the  fields  at  the  boundary  of  the  (m-l)th  and  mth 
media  are  given  by 


where  the  superscripts  +  and  -  designate  waves  traveling  in  the 

forward  and  backward  z  directions.  r  and  t  are  the  Fresnel 

m  m 

coefficients  and  b  is  the  optical  phase,  given  by  6  =  2nn  d  /X 

m  m  m  m  o 

where  dm  is  the  layer  thickness.  This  matrix  equation  is  a  recur¬ 
sion  relation  between  the  electric  fields  in  successive  layers  and 
allows  convenient  computation  of  E(z)  in  terms  of  the  incident 

field  E+  . 
o 

The  indices  of  refraction  for  the  film  and  substrate 
materials  at  the  ruby  laser  wavelength  that  were  used  in  our  calcula¬ 
tions  are  listed  in  Table  5. 

Solutions  were  obtained  for  normal  incidence  on  single¬ 
layer  films,  bilayer  V-type  antireflection  coatings  and  multilayer 
reflection  coatings.  Solutions  for  non-normal  incidence  on  single¬ 
layer  coatings  for  polarizations  perpendicular  (S)  and  parallel  (P) 
to  the  plane  of  incidence  were  also  obtained.  A  brief  summary  on 
general  expressions  of  (  E./  E^  |  for  different  cases  as  well  as  the 
corresponding  computer  programs  can  be  found  in  Appendix  A.  Some 
of  the  results  from  the  calculations  will  be  used  to  analyze  our  experi¬ 
mental  data  in  the  subsequent  sections. 
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Table  5.  Refractive  indices  of  film  and  substrate  materials  at 
6943A. 


Material 

Index  of  refraction 

Reference 

MgF? 

1.  38 

[9] 

sio2 

1.456 

[9] 

ThF . 

4 

1.  52 

[10J 

Zr02 

1.975 

[93 

Ce02 

2.  2 

[n3 

Ti°2 

2.28 

[93 

ZnS 

2.32 

[103 

Fused  silica 

1.456 

[123 

BSC2  glass 

1.  513 

[133 

NaCl 

1.  54 

[l  4J 

Spinel 

1.73 

[153 

YAG 

1.83 

[153 

a.  Single -layer  Films 

In  the  case  of  normal  incidence.  Figure  23  shows  the 
distributions  of  the  electric-field  intensity  for  quarter-wave  films  of 
MgF^»  SiO^*  ZrC>2  and  TiO^  on  a  glass  substrate.  Note  that  the 
intensity  scale  is  broken  at  *.6,  and  the  peak  of  some  distributions 
in  air  has  been  displaced  downward  and  plotted  on  the  scale  in 
parenthesis.  Figure  24  shows  the  intensity  distributions  for  half¬ 
wave  films  of  the  same  materials.  Since  Si02  films  have  an  index 
of  refraction  near  to  the  index  of  the  glass  substrate,  there  is  not 
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much  difference  between  the  intensity  distributions  for  the  quarter 
and  half-wave  films.  The  main  differences  for  the  other  film  mater¬ 
ials  occur  at  the  air-film  interface. 

In  the  case  of  non-normal  incidence,  Figures  25  and  26 
show  the  relative  power-density  distributions  for  both  S  and  P  polari¬ 
zations  for  quarter  and  half-wave  thick  films  of  MgF The  thres¬ 
holds  for  P  polarization  should  be  significantly  lower  than  those  for 
S  polarization. 

b.  Bilayer  Antireflection  Coatings  (V-type) 

The  electric  fields  in  bilayer  antireflection  coatings  of 
a  V-design  were  calculate:  d  for  normal  incidence  in  two  steps.  First, 
the  thicknesses  of  the  two  layers  were  calculated,  and  then  the  field 3 
were  computed.  The  thicknesses  were  calculated  directly  using  the 
relations  of  Hass  [16]  and  Catalan  [17].  These  relations  yield  the 
thicknesses  required  to  obtain  0%  reflectance  at  a  given  wavelength 
for  normal  incidence.  The  actual  coatings  were  designed  for  pairs 
of  films  with  high  and  low  refractive  indices  on  BSC-2  glass  substrates. 
The  higher-index  layer  was  adjacent  to  the  substrate  and  was  less 
than  X/4  thick.  The  calculated  optical  phases  and  thicknesses  and 
the  physical  thicknesses  are  listed  in  Table  6.  The  substrate  was 
BSC-2  glass.  The  thicknesses  of  the  high-index  layers  ranged  from 
1/5  to  1/2  of  a  quarter-wave,  and  the  thicknesses  of  the  low  index 
layers  were  all  about  1.3  quarter -waves. 

The  electric-field  distributions  shown  in  Figure  27  are 
about  the  same  for  each  V-coat.  The  only  differences  occur  in  the 
power  density  distributions  which  is  obtained  by  multiplying  the 
relative  electric-field  intensity  by  the  index  of  refraction. 


Index 

Substrate 

Film 

BSC-2Glass 

1- Zr02/Si02 

2- Ti02/MgF2 

3  Zr02/MgF2 

4  Ti02/Si02 


incidence 


Figure  27.  Relative  electric-field  intensity  and  power-density 
distributions  for  bilayer  antireflection  coatings  or. 
BSC-2  glass  substrates. 


Table  6.  Thicknesses  of  bilayer  antireflection  coatings. 


Sample  No.  0143  0145  0147  0142 


Materials 
6^  (degrees) 
&2(degrees) 
n j d j  ( X /  4) 
n2d2a/4, 

djfnm) 

d2(nm) 


Zr02/Si02 
111.9 
40.6 
1.  244 
0.451 
148.  2 
39.65 


ZrOz/MgF2 
112.  5 
29.  55 
1.  250 
0.  328 
157.  2 
28.86 


TiOz/MgF2 
116.65 
19.  2 
1.296 
0.  214 
163.  0 
16.  27 


Ti02/Si02 
118.8 
25.2 
1.320 
0.  280 
157.3 
21. 29 


c.  Multilayer  reflectors 

The  electric  field  distributions  for  normal  incidence  in 
two  multilayer  reflection  coatings  arc  shown  in  Figures  28  and  29, 
and  the  fields  at  the  film  interfaces  of  these  and  other  reflectors 
are  listed  in  Table  7.  Note  that  the  most  damage-resistant  reflectors 
of  the  Ti02/Si02  configuration  have  lower  electric  fields  at  the 
H-L  interfaces  than  the  other  designs.  Furthermore,  for  all  reflectors 
the  fields  decrease  rapidly  from  the  top  layer  inward,  and  at  the 
substrate  they  are  negligible.  This  explains  why  only  the  top  layers 
are  damaged  at  threshold. 

As  illustrated  in  Figures  28  and  29  the  electric  field 
intensities  peak  at  the  H-L  interfaces.  Therefore,  for  laser  irradia¬ 
tion  above  threshold,  damage  involves  removal  of  layers  in  pairs 
since  the  high-index  films  are  less  damage  resistant  than  the  low- 
index  films.  For  example,  at  the  first  peak  in  the  electric-field- 
squared  inside  the  Ti02/SiC>2  reflector,  the  energy  density  in  the 
high  index  layer  was  65%  higher  than  in  the  low-index  layer. 


Incidence 


Figure  29. 


Relative  electric-field  intensity  distributions  for  a 

10  2 

22-layer  TiC>2/Si02  reflector  with  G(HL)  HL  A 
configuration. 


Table  7.  Normallaed  electric  field  Inleneltlee  In  multlleyer  r. Hector. 


Sample 

Film 

Material* 

Dealgn 

reflectance 

Number 
of layera 

Film  order 

lE'<l ' 

In  top 
low  Index 
layer 

i*/e;i2 

at  (Irat  II-  L 
Interface 

|E/<|2 
at  eecond  II.  L 

% 

0102.  10) 

TI02/Si02 

99.  9* 

20 

G(HL)9IIL*A 

1.  887 

0.724 

0.278 

OIOS  10b 

Zr02/SIC2 

99.9* 

29 

C(IIL)UHA 

-- 

I.02S 

0.  SS7 

0107.  108 

Zr02/M*F2 

99.  9* 

2S 

G(ISL)I2MA 

-- 

I.02S 

0.  SOI 

o  o 

M  i— 

o 

IS* 

c«o2/sio2 

9S.  ) 

1) 

ohu6iia 

-- 

1.027 

0.  S7) 

OIO).  104 

Ti02/SI02 

99.9 

22 

C(HU,0IIL2A 

1.  887 

0.769 

0)14 

0109.  110 

98.7* 

12 

cihd'iil2a 

1. 87S 

0.764 

0.  )  1 2 

None 

98.7 

II 

C(IIL)*!IA 

•• 

0.764 

0.  )  1 2 

None 

•  1 

97.  ) 

12 

G(IIL|4A 

).  94S 

1.609 

0.6S6 

Incidence  from  eubetrate  aide: 

OIO).  104 

tio2/sio2 

Aa  above 

-- 

-- 

1.  116 

0.4SS 

OIOS.  106 

Zr02/SI02 

Aa  above 

-- 

-- 

1.487 

0.  808 

CI07.  108 

Zr02/M*F2 

Ae  above 

•• 

I.48S 

0.72S 

3.  Damage  thresholds  for  entrance  and  exit  face  coatings. 

The  damage  thresholds  of  entrance  and  exit-face  coatings 
were  measured.  These  tests  were  conducted  to  find  out  if  there 
is  a  direct  relationship  between  the  thresholds  and  the  electric  field 
standing -wave  distributions  in  the  films.  Listed  in  Table  8  are  the 
peak  energy  density  thresholds  as  measured  by  LIS  for  the  s  ingle - 
layer  films  and  by  spark  detection  for  the  multi-layer  reflectors. 

In  every  case,  except  for  the  ZrC>2/MgF2  reflector  (sample  0108), 
the  entrance-face  thresholds  were  equal  to  or  greater  than  the  exit- 
face  thresholds. 


T»l>lr  R.  Damage  tliresluilds  (or  entrance  and  exit  fair  coaling* 


Sa  tuple 

Coating 

description 

Spot-  'me 

l»eak  energy  density* 
entrance  exit 

Ratio 

lent  r.  /exit ) 

Single  layers 

nun 

J  /cm^ 

J /cm* 

0125 

Si02. 1/4 

0.  055 

1  17 

1 10.-.II4 

1.0 

OIJ1 

MgF2.»/4 

0.  055 

III. Ml 

•>124 

1 . 0 

Ol  J5 

MgF2.)/2 

0.07 

65. 102 

44 

1.9 

0121 

TIOr  »/4 

0.  072lrntr.  ) 

0.  09  (exit) 

)6.  42 

20-21 

1.9 

01)1 

ZrOj  Ml 

0.  072(cnt  r.  ) 

0.  06  Jlcxit) 

12-1) 

9.5-10 

1.  ) 

01)2 

ZrO  l.\/l 

0.  15 

a.  5-9 

8 

1.  1 

Multilayer 

reflectors 

OIO) 

TiC2/SI02 

0.  056(cntr.  ) 
0.07  (exit) 

121-126 

95-116 

1.  2 

0 104 

tio2/sio2 

0.056 

117. 127 

98- 108 

1.2 

0106 

ZrO2/SI02 

0.  1) 

18.  5 

14.  5 

1.  ) 

0108 

ZrC2/MgF2 

0.  062(entr. ) 

0.  055(exit) 

90 

128-1)1 

0.7 

“Single  layrr  threshold*  were  determined  by  laser- Induced  scattering  (LIS).  Reflector 
thresholds  were  determined  by  spark  formation. 


For  these  tests,  each  coating  was  alternately  faced 
toward  and  away  from  the  incident  laser  beam.  This  result  is 
similar  to  that  generally  observed  for  glass  and  crystal  surfaces  [l8J, 
but  there  has  been  considerable  discussion  about  which  surface, 
front  or  back,  at  which  a  laser-induced  optical  discharge  first  occurs. 
Fersman  and  Khazov  [19J  and  Dupont  et  al.  [20]  in  measuring  surface 
damage  of  glass,  observed  a  spark  first  at  the  front  surface,  although 
damage  effects  at  the  rear  surface  were  more  severe.  In  further 
investigations  Giuliano  [2 1 3  detected  microscopic  damage  at  the  exit 
surface  of  sapphire  prior  to  the  development  of  a  visible  spark. 

This  latter  result  is  in  agreement  with  the  present  observations  for 
thin  film  coatings  in  which  damage  detectable  by  laser-induced 
scattering  (LIS)  occurred  before  or  at  the  spark  threshold.  Now  since 
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the  single-layer  thresholds  were  determined  by  LIS,  any  possible  con¬ 
fusion  in  damage  detection  by  sparks  is  avoided.  However,  there 
might  be  a  difficulty  in  the  comparison  of  multilayer  thresholds  since 
their  thresholds  were  determined  by  spark  formation. 

The  electric-field  distributions  for  exit-face  single¬ 
layer  films  of  quarter  and  half-wave  thicknesses  are  shown  in  Figures 
^0  and  31.  These  distributions  can  be  compared  to  those  for  entrance- 
face  coatings  presented  in  Figures  23  and  24.  The  obvious  difference 
is  that  maxima  invariably  occur  at  the  air-film  interface  for  each 
exit-surface  coating.  When  the  quarter-wave  coating  is  moved  from 
the  entrance  to  the  exit  face,  the  intensity  increases  at  the  air-film 
interface  and  decreases  at  the  film-glass  interface  with  the  total 
intensity  within  the  film  slightly  increasing.  For  TiO  ,  the  intensity 
increases  by  3.  55  at  the  air-film  interface  and  decreases  by  0.  30 
at  the  fL/n-glass  interface.  The  measured  ratio  of  entrance  to  exit- 
face  damage  thresholds  was  1.9  for  TiO^  this  value  corresponds 
well  to  the  calculated  ratio  of  1.  56  for  the  peak  intensities  (at  the 
film-glass  interface  in  the  entrance  orientation,  and  at  the  air-glass 
interface  in  the  exit  orientation).  For  MgF^  and  SiOz,  the  measure^ 
ratio  of  entrance-to-exit  thresholds  was  unity,  corresponding  to  the 
calculated  ratio  of  the  total  intensities  within  the  film  (1.09  and  1.04, 
resp. ). 

When  the  half-wave  films  are  considered,  there  is  a 
simplification  in  that  the  film  material  does  not  affect  the  entrance- 
to-exit  intensity  ratios.  When  moved  from  entrance-tc-exit  face,  the 
intensity  at  the  air-film  interface  increases  by  1.58  regardless  of  film 
material.  The  intensity  at  the  film-glass  interface  is  the  same  as  the 
air-film  interface,  while  the  intensity  extremum  within  the  film 
decreases  by  0.69.  For  MgF2  and  ZrC>2  half-wave  films,  the 
measured  ratios  of  entrance-to-exit  thresholds  were  1.9  and  1.3. 

The  average  of  these  measured  ratios  is  1.6,  quite  close  to  the  value 
for  the  air-film  interface. 


E  Substrate  Film 

ir  Fq  BSC-2Glass  ^  Thick  Air 


Figure  30.  Relative  electric-field  intensity  distributions 
for  quarter-wave  films  coated  on  the  exit 
face  of  BSC-2  glass  substrates. 


Incidence 


Figure  31 


Relative  electric-field  intensity  distributions  for 
half-wave  films  coated  on  the  exit  face  of  BSC -2 
glass  substrates. 
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The  electric-field  intensity  distribution  in  an  exit- 
surface  TiO^/SiO^  reflector  is  shown  in  Figure  32  and  may  be  com¬ 
pared  to  the  distribution  in  Figure  29  for  the  same  reflector  in  the 
usual  front-surface  orientation.  The  maxima  and  minima  occur  at 
the  layer  interfaces  for  both  orientations,  but  higher  peaks  are  pre¬ 
sent  in  the  reversed  reflector.  Another  difference  is  that  the  extreme¬ 
ly  high  peaks  of  the  standing  waves  which  occur  in  front  of  the  reflec¬ 
tor  in  Figure  29  are  relocated  in  the  glass  substrate  in  Figure  32. 

The  peak  values  of  the  intensity  are  45%  higher  for 
the  reverse  orientation  and  therefore  by  elective  field  considerations 
-lone,  the  damage  thresholds  would  be  expected  to  be  45%  less. 
However,  the  spark  thresholds  of  reversed  TiO  /SiO  and  ZrO  /SiO 

m  2  2  2 

reflectors  were  only  20%  and  30%  lower  and  the  ZrO^/  MgF^ 
reflector  was  about  45%  higher!  Apparently  the  substrata  caused 
no  degradation  of  the  threshold.  Quite  to  the  contrary,  the  sub¬ 
strate  probably  provided  structural  support  to  the  adjacent  layers 
which  experienced  the  highest  electric  fields.  Furthermore,  the 
residual  stresses  in  the  layers  near  the  substrate  theoretically  have 
less  residual  tensile  stress  than  the  layers  near  the  air-film  inter¬ 
face.  It  is  then  reasonable  that  the  layers  near  the  substrate  are 
more  damage  resistant  for  a  given  energy  density. 


Tf  ( 
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V.  ROLE  OF  COATING  DEFECTS  IN  LASER -INDUCED 
DAMAGE  TO  THIN  FILMS 

1.  Introduction 

It  is  noted  that  in  the  literature  the  spot-size  of  the 
laser  beam  was  seldom  reported  along  with  its  corresponding  laser- 
induced  damage  thresholds  of  laser  materials.  One  of  the  most  impor¬ 
tant  experimental  results  in  correlating  laser  parameters  in  this  study 
is  that  the  damage  threshold  of  the  thin  film  increases  as  the  spot- 
size  of  the  laser  beam  decreases.  An  80%  variation  of  damage  thres¬ 
hold  for  the  same  film  could  be  observed  by  just  varying  the  spot-size. 

In  this  report  we  develop  a  simple  model  correlating 
the  nature  and  distribution  of  coating  defects  to  this  spot-size  depen¬ 
dence.  The  model  assumes  that  damage  due  to  local  defects  or  impuri¬ 
ties  in  the  film  is  different  from  the  intrinsic  damage  of  the  material 
and  that  the  defect  sites  are  randomly  distributed  on  the  coating  surface. 
The  probability  of  the  laser  beam  striking  a  defect  site  with  a  certain 
intensity  can  be  calculated  when  the  intensity  distribution  of  the  laser 
beam  is  known. 


2.  The  Defect  Model 

For  a  random  distribution  of  points  on  a  plane,  the  pro¬ 
bability  that  a  randomly  chosen  area  (e.  g. ,  a  circle  of  radius  r)  will 
contain  exactly  n  points  can  be  described  by  the  Poisson  function 

P(n)  =  (  (  pnr2)n  /  n !  }  exp  (-pnr2)  (5) 

where  p  is  the  mean  surface  density  of  the  defects.  The  probability 
that  no  points  will  be  contained  in  this  area  is 

P(0)  =  exp  (  -  pnr2 )  (6) 
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which  is  the  probability  that  the  area  will  contain  no  point  within  a 
distance  r.  Hence,  the  proportion  of  distances  of  nearest  neighbors 
less  than  or  equal  to  the  distance  r  is 

p(r)  =  I  -  exp(-pnr2)  .  (7) 

Now,  if  r  is  allowed  to  vary,  the  probability  distribution  of  r  is  then 
expressed  as 

dP(r)  =  2prrr  exp  (  -  pnr“)  dr  .  (8) 

Thus,  if  we  consider  the  case  of  a  square  pulse  with  width  w  less 
than  or  equal  to  r,  the  probability  of  such  a  pulse  hitting  one  or  more 
defects  will  be  easily  derived  as 

P(w)  =  1  -  exp [  -  (w/d  )2]  (9) 

“  o 

where  dQ  is  the  expectation  value  of  r,  or  the  mean  distance  of  two 

defects,  and  is  related  to  the  mean  density  p.  In  the  actual  case, 

a  TEM00  ruby  laser  pulse  was  used  so  that  the  square  pulse  should 

be  modified  to  the  form  of  a  Gaussian  beam  with  spot-size  w  .  Then 

o 

Eq.  (9)  becomes 

P(wq)  =  1  -  exp [  ~  2  (wQ/do)2J  (10) 

Now,  by  assuming  that  the  coating  defects  on  a  dielec¬ 
tric  filrr  surface  have  a  damage  threshold  1^  (in  joules/cm2),  much 
less  than  the  intrinsic  damage  threshold  Ij  of  the  film,  the  total 
damage  threshold  I  can  be  expressed  in  terms  of  the  probability  given 
in  Eq.  (10). 


I 


IdP,wo)  +  '|  l  1  -  P(wo>  J 
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To  analyze  our  experimental  data,  it  is  convenient  to  normalize  the 
total  damage  threshold  I  with  respect  to  the  defect  damage  threshold 
such  that 


J  =  1  +  (T|  -  1 )  exp [  -  -  a*1  2  (w  /d  )2]  (11) 

o  o  o 

where  we  have  defined  the  ratios  J)=I/I  and  Tj  =  I .  / 1  .  If  the  distribu- 

u  in 

tion  of  defects  on  the  coating  is  random,  only  one  parameter,  namely 
d  ,  is  needed  to  completely  describe  the  distribution.  1^  and  T}  will 
be  parameters  which  can  be  related  to  the  type  of  defect,  whether 
defect  damage  is  associated  to  absorption  or  electron  avalanche. 

In  a  digression,  let  us  take  a  closer  look  at  the  para¬ 
meter  T] ,  which  can  be  viewed  as  a  function  of  the  defect  damage 
mechanism.  A  recent  calculation  by  Bloembergen  [22]  concerning  the 
role  of  pores  in  laser-induced  damage  to  surfaces  can  be  extended  to 
structural  defects  in  films.  If  a  defect  was  a  cylindrical  groove  on 
the  surface,  then 

=  4n4/  (n^  +  1)^  (12) 

which  has  a  maximum  value  of  4  when  the  refractive  index  is  infinite. 
If  a  defect  was  a  needle  cavity  or  the  so-called  'crack'  defect,  which 
is  an  oblate  ellipsoidal  void  in  the  surface,  the  enhancement  factor  Tl 
will  have  a  much  larger  value  of 

T|  ~  n4.  (13) 

Other  types  of  defects  will  have  different  values  of  the  enhancement 
factor.  Hence,  the  parameter  r\  is  very  sensitive  to  the  type  of  defect 
on  the  coatings  or  the  surface  of  the  bulk  materials.  Table  9  gives 
some  estimates  of  the  enhancement  factor  B  for  five  different  films. 

To  illustrate  the  results  of  such  a  probability  model, 
Figure  33  displays  a  theoretical  plot  of  the  total  damage  threshold  I 
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Figure  33.  Theoretical  plot  of  the  total  damage  threshold 
versus  spot-size. 


■ 


I 
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versus  spot-sizes  wq.  For  very  large  spot-sizee,  the  measured 

threshold  will  be  the  defect  damage  threshold  1^.  For  decreasing  spot- 

sizes.  the  damage  threshold  will  follow  a  variety  of  curves  depending 

on  the  distribution  of  defects,  e.  g.  ,  the  mean  distance  d  .  In  the  case 

o 

of  thin-film  coatings,  the  damage  thresholds  are  very  low  because  d 

o 

is  a  small  fraction  of  the  spot-sizes  used  in  the  usual  laser  damage 

tests.  Hence,  the  curve  i6  similar  to  the  one  for  d  =  0.  5  w  .  and 

o  o 

the  threshold  value  is  close  to  the  defect  damage  I,  and  not  to  the 

d 

intrinsic  damage,  I  .  However  in  the  case  of  single  crystal  damage 
studies,  e.  g.  ,  for  NaCl  [23]  the  spot-size  is  very  small  with  respect 
to  the  defect  distance  d^.  Then  the  curve  is  similar  to  the  one  for 
V  4wo'  and  in  fact*  the  measured  threshold  would  be  very  close  to 
the  intrinsic  damage  I  .  Therefore,  using  this  model  both  thin  film 
and  surface  damage  can  be  related  for  a  single  material. 

Table  9.  Some  estimates  of  enhancement  factor  r|  due  to  structural 
defects 


Sample 

n 

Groove 

T|  =  4n4/  (n2  +  1  )2 

Crack 
T|~  n4 

ZnS 

2.  32 

2.  84 

29 

TiOz 

2.  28 

2.81 

27 

ZrOz 

1. 98 

2.  54 

15 

sio2 

1.46 

CO 

w 

4.  5 
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3.  Experimental  Results 

The  descriptions  of  the  experimentai  setup  (see  Fig.  3) 

and  procedure  used  in  this  experiment  can  be  found  in  Section  II  and 

III.  Three  methods  were  used  to  monitor  the  onset  of  the  damage: 

spark  detection,  optical  microscopy  and  laser-induced  scatter  (LIS). 

However  here  we  emphasize  that  the  LIS  technique  was  our  primary 

tooi  for  diagnosis.  Spot-sizes  ranging  from  52(am  to  280|im  were 

used  and  areas  of  a  test  sample  were  irradiated  by  one  shot  only. 

The  mean  distance  d^  between  defects  of  each  sample  was  determined 

by  measuring  the  density  of  coating  defects  for  ten  randomly  chosen 

areas  on  the  sample,  which  was  examined  under  a  scanning  electron 

micro  cope  with  magnification  of  1000. 

A  half-wave  film  of  ZrO  and  a  quarter-wave  film  of 

£ 

ZnS  with  glass  substrate  were  used  in  studying  the  spot-size  depen¬ 
dence  on  damage  thresholds.  By  varying  the  spot-size  of  the  laser 
beam,  damage  thresholds  were  determined  by  the  onset  of  LIS.  The 
laser  pulsewidth  was  8  nsec.  Figures  34  and  35  display  the  experi¬ 
mental  data  as  well  as  the  theoretical  curves  using  the  forementioned 

defect  model.  For  ZrO  film,  the  defect  damage  threshold  I.  was 

c  2  d 

found  to  be  9.  2  joules/cm  and  the  data  was  fitted  to  a  curve  with 

ri  of  3  and  dQ  of  20,4  microns.  Figure  36  compares  damage  data 

for  three  other  coatings  to  the  curves  derived  from  the  defect  model. 

For  the  plotting  of  this  data,  dQ  was  measured  by  SEM  photographs 

and  was  not  fitted.  The  sketchy  nature  of  the  data  for  these  three 

coatings  prevented  an  accurate  determination  of  the  enhancement 

factor  T| ,  but  it  is  to  be  noted  that  the  data  was  consistent  with  the 

model. 


Spot-size  dependence  for  single  half-wave  thick  film  of  ZrO £  on  BSC- 
glass  substrate.  j?  =  I/I,  and  ti  =  I. /I,. 


•  MgF2  ,X/4(d0=47/xm,Id  =  l05J/cm2) 

A  Zr02  ,X/2(d0=  30/xm,Id  =  l5.5  J/cm2) 

■  Ti02/Si02,22!oyer(d0=25Atm,Id=73J/cm2) 


SPOT- SIZE  WQ  ( in  units  of  d0) 

Figure  36.  Normalized  damage  threshold  as  a  function  of 
laser  beam  spot-size  including  the  role  of 
local  defects.  The  range  of  opofc-size  is  55 
to  200  microns. 
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4,  The  Roie  of  Coating  Defects 

A  closer  look  at  the  morphology  of  the  coating  defects 

or  impurities  on  the  films  using  the  electron  microscope  has  indicated 

that  defects  are  irregular  both  in  size',  :nd  shapes  and  the  density 

may  strongly  depend  on  several  film  parameters  and  preparation 

methods.  In  general,  six  types  of  defects  were  observed.  Figures 

37  to  43  show  that  such  local  defects  can  be  found  in  different  films. 

For  the  ZnS  film  used  in  the  experiment,  five  types  of  defects  were 

found  (see  Fig.  44).  It  is  recalled  that  the  measured  d  was  20.4 

o 

microns  and  the  fitted  d  was  30  microns,  which  indicates  that  the 

o 

damage  may  be  primarily  caused  by  those  defects  having  irregular 
forms  (or  "mountain  chains")  and  "shadows  with  strips".  One  can 
also  easily  observe  in  this  case  that  a  small  mean  distance  between 
defects  occurs  for  the  most  prevalent  defect,  "shadows  with  a  dot". 

By  examining  coating  procedures,  one  may  be  able  to 
change  the  mean  distance  between  coating  defects.  At  USC,  we  have 
l 'tempted  to  vary  the  coating  defects  on  one  type  of  film  by  varying 
the  deposition  rate  in  the  preparation  of  the  coatings.  It  was  fotud 
that  by  doing  so,  dQ  decreases  as  the  deposition  rate  increases. 

Figure  45  simply  illustrates  this  result  for  MgF  films.  The  idea 

Ct 

was  then  to  damage  all  three  of  these  films  of  the  same  material 
and  find  out  how  it  relates  to  the  parameter  d  . 
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Figure  37.  Scanning  Electron  Microscope  (SEM)  photograph  of 
defects  of  the  two  extreme  forms,  a  hill  (upper)  and 
a  hole  (lower),  on  an  eleven  layer  CeO^/SiO^  reflector. 
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Figure  38.  SEM  photograph  of  defcctB  of  a  group  of  holes  on  a 
single  quarter-wave  MgF^  film. 


Figure  39.  SEM  photograph  of  defects  in  the  form  of  a  crack  on 
single  quarter-wave  MgF^  film. 


SEM  photograph  of  defects  in  the  form  of  irregular 
shaped  hills,  or  "mountain  chain,  "  on  a  single 
quarter-wave  MgF2  film. 


SEM  photograph  of  defects  in  the  forms  of  shadow 
with  dot,  crack  and  hill  on  a  single  quarter-wave 
ZnS  film. 
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«  44.  A  schematic  plot  of  size  versus  distribution 
and  nature  of  coating  defects,  for  a  single 
quarter-wave  film  of  ZnS  (sample  number  U5a). 
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VI.  OTHER  LASER  PARAMETERS 


As  described  in  Section  II,  knowledge  of  the  spatial  and 
temporal  behavior  of  the  laser  beam,  and  hence  the  mode  properties 
of  the  laser,  is  required  for  an  accurate  prediction  of  damage. 
Furthermore,  if  we  want  to  know  whether  or  not  time -dependent 
damage  mechanisms  are  involved,  correlation  of  the  damage  threshold 
to  the  laser  pulsewidth  should  be  studied.  In  this  program,  two  experi¬ 
ments  were  performed  to  investigate  thin-film  damage  as  a  function 
of  these  two  laser  parameters,  namely,  the  transverse  mode  struc¬ 
ture  and  the  pulsewidth  of  the  laser. 

1.  Damage  by  Multimode  Lasers 

For  the  mu.'timode  damage  Btudy,  the  ruby  laser 
oscillator  was  adjusted  to  have  a  Fresnel  number  of  8.  0  instead  of 
0.4,  which  allowed  the  oscillator  to  operate  in  many  transverse 
modes.  For  our  laser,  the  cavity  FrcBnel  number  of  8.0  corresponded 
to  an  aperture  having  a  diameter  of  4.  47  mm.  The  test  samples  were 
placed  at  approximately  20  cm  after  a  31 -cm  focal  length  lens  located 
at  116  cm  from  the  aperture.  In  single-mode  operation,  the  lens 
would  see  a  far-field  diffraction  pattern  at  this  distance  (1.3R)  [24] 
but  it  was  in  the  near-field  distances  of  the  samples  from  the  aperture 
were  computed  to  be  approximately  0.1 25R  by  using  the  geometrical 
lens  transformation. 

The  damage  thresholds  (spatial  averages)  of  these 
dielectric  reflectors  (samples  S101,  102  and  103)  are  listed  in  Table 
10  for  near-field,  non-single  transverse-mode  laser  irradiation. 

Shown  for  comparison  are  the  average,  single-mode  thresholds  com¬ 
puted  from  the  data  given  in  Table  4  of  Quarterly  Technical  Report 

No.  4.  An  "average  energy-density  threshold  for  a  Gaussian  trans- 

2 

verse  intensity  is  defined  by  E=E^,/nw  .  This  average  is  equal  to 
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one-half  the  axial  energy  density  prescribed.  Because  multimode 
emission  is  temporally  unstable,  theoretical  computations  by 
McAllister  [25]  of  the  time-dependent  intensity  distributions  plus 
geometrical  lens  optics  were  used  to  estimate  the  spot-size  at  the 
sample  locations.  Since  these  measurements  were  performed  prior 
to  incorporating  the  spark  detection  system,  the  tl  resholds  of  the 
reflectors  were  established  using  low-power  (7  X  and  40  X  )  microscopes. 
Taking  into  account  the  differing  sensitivity  of  the  detection  methods, 
the  multimode  thresholds  were  still  considerably  less  than  those  for 
single-mode  pulses.  A  typical  comparison  of  damage  morphology 
due  to  these  two  different  pulses  can  be  found  in  Figure  46. 


Table  10.  Damage  thresholds  for  non-single  transverse-mode  laser 
irradiation  (nominal  pulsewidth  is  16-18  nsec) 


Sample 

Description 

Spot-size, 

Average 

Average  energy 
density 

Peak  power 
Density,  avg. 

mm 

J/cm^ 

GW  /cm" 

S101 

Ti02/Si02 

0.  65 

10  (22-28) 

0.  6(1.  5-2.  0)a 

SI  02 

Ti02/Si02 

0.  8 

7  (25) 

0.4(1. 8) 

SI  03 

Zr02/Si02 

1.2 

3  (4) 

0.  14(0.  26) 

For  comparison,  damage  thresholds  for  single-mode  laser  irradiation 
are  given  in  parentheses. 


The  damage  thresholds  of  thin  film  coatings  for  multi-trans- 
verse-mode  (MM)  laser  radiation  were  much  lower  than  for  singlc-trans - 
verse-mode  (SM)  radiatior.  For  the  three  multilayer  samples  listed  in 
Table  10,  the  energy  density  thresholds  for  MM  radiation  were  40%, 

30%  and  75%  of  the  respective  SM  thresholds.  Apparently,  the  more 
damage  resistant  the  coating,  the  greater  the  difference  between  MM 


Figure  46.  Laser-induced  damage  far  above  threshold  on  a 
TiC^/SiOg  multilayer  reflector  (Sample  S102) 
caused  by  the  output  of  (a)  a  single-mode  oscillator 


and  (b)  a  muLl-transverse-mode  oscillator.  Damage 


at  film  defects  is  apparent.  Photographed  with  dark 


field  Illumination.  Magnification:  1  20  X  . 
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and  SM  thresholds.  It  is  only  apparent  that  the  ratios  for  MM  to  SM 
powe r  density  thresholds  are  even  lower  than  those  for  energy  density. 
When  the  longer  pulsewidths  used  in  the  MM  tests  are  taken  into 
account,  the  ratios  are  about  the  same  for  both.  Theoretical  calcula¬ 
tions  of  the  MM  output  have  shown  that  spatial  nonuniformities  in  the 
power-density  distribution  at  a  given  instant  in  time  can  be  consider¬ 
ably  more  severe  than  for  the  energy-density  distribution  which  is  an 
integration  over  the  pulse  duration.  A  larger  difference  for  power- 
density  thresholds  would  have  indicated  that  the  damage  mechanisms 
are  strongly  time  dependent  in  this  pulsewidth  range.  Thus  the 
experimental  thresholds  indicate  the  contrary. 

The  calculation  of  the  MM  thresholds  was  very  approxi¬ 
mate  because  the  spatial  distribution  of  the  output  at  the  sample  plane 
was  both  radially  nonuniform  and  time  dependent.  Numerical  calcu¬ 
lations  of  the  theoretical  intensity  distributions  at  three  different 
times  during  a  single  pulse  were  supplied  by  McAllister  [25]  and  are 
shown  in  Figure  47.  Each  curve  is  normalized  to  its  peak  radial 
value,  and  radial  distances  are  normalized  to  the  aperture  radius  a. 
The  curves  represent  the  output  in  the  plane  of  the  irradiated  samples 
which  was  C.  125R  from  the  laser.  This  distance  is  in  the  near  dif¬ 
fraction  field  at  which  the  oscillator  aperture  subtended  four  Fresnel 
zones  as  viewed  by  the  samples. 

The  temporal  and  spatial  structure  of  this  MM  output  was 
semi-controlled  by  use  of  an  aperture  which  set  the  oscillator  Fresnel 
number  N  equal  to  8,0.  Without  the  aperture,  N  would  have  equaled 
16  as  established  by  the  diameter  of  the  ruby  rod,  and  an  even  more 
erratic  temporal  and  spatial  output  would  have  resulted.  Additional 
computer  analysis  of  the  laser  output  at  several  distances  in  the  neur 
and  far  fields  revealed  that  the  spatial  nonuniformities  at  a  given 
instant  decrease  with  increasing  dietance,  but  temporal  variations 
persist.  For  example,  at  one  Rayleigh  distance  '7.2  meters  from  the 
4.47mm  aperture)  the  radial  distribution  is  Gaussian-like  during  the 
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Flgvjre  47.  Theoretical  transverse  intensity  distributions 

for  the  output  of  a  multi-transver sc-modc  ruby 
laser  (N=8.  0)  at  0.  125  R  from  the  oscillator 
aperture.  Curve  1  (3)  occurs  during  the  rise 
(fall)  of  the  temporai  pulse  when  the  total 
power  equals  0.  51  (0.  55)  the  peak  power. 

Curve  2  occurs  at  the  pulse  peak. 


L 


entire  pulse,  but  the  spot-«ize  severely  decreases  in  time  after  the 
pulse  peak. 

The  damage  thresholds  of  thin  film  coatings  irradiated 

by  MM  later  pulses  were  considerably  lower  than  with  SM  pulses. 

This  is  in  agreement  with  the  results  of  two  other  investigators  for 

other  dielectric  surfaces.  lzamitani,  Hosaki  and  Yamanaka  [26] 

measured  the  single-shot  MM  surface  damage  threshold  of  Hoya  glass 

to  be  about  one-half  the  value  of  the  SM  threshold  which  was  28  J/cm2. 

Similarly,  Bass  [27]  determined  that  the  surface  damage  threshold  of 

LiNbO  caused  by  a  single  pulse  from  a  multimode  NdrYAG  laser 
2 

was  0.8  GW/cm  (average)  in  contrast  to  that  caused  by  a  SM  pulse 

2 

of  average  power  density  equal  to  1.6  GW/cm  .  The  lower  damage 
resistance  to  MM  pulses  is  in  agreement  with  the  general  results  of 
our  thin  film  experiments.  Relative  "hot  spots"  in  the  MM  transverse 
beam  profile  or  excitation  of  additional  damage  mechanisms  by  en¬ 
hanced  temporal  and  spatial  gradients  may  cause  early  damage.  It  is 
not  clear  that  these  arc  the  only  possible  explanations.  The  difficulty 
in  performing  accurate  calculations  of  the  MM  power  may  well  be 
responsible  for  some  of  the  apparent  difference  [28]. 

2.  Damage  by  Different  Pulsewidth  of  the  Laser  Beam 

Single,  half-wave  layer  ZrO.,  films  were  irradiated  by 

the  ruby  laser  using  11  to  35  nsec  pulsewidths.  The  damage  thresholds, 

shown  in  Figure  48,  exhibited  an  interesting  dependence  on  pulse 

duration.  The  threshold  detected  by  LIS  increased  by  36%,  from  12.  5 

2 

to  17  J/cm  over  the  pulsewidth  range  of  13  to  32  nsec.  Over  the 
same  temporal  range,  the  spark  threshold  rose  from  17.5  to  30  J/cm2, 
an  increaoe  of  67%.  The  slopes  of  the  straight  lines  drawn  between 
the  data  points  lie  between  the  hypothetical,  constant  energy-density 
curves  (horizontal  lines)  and  the  constant  power-density  curves  Grawn 
through  the  shorter  pulsewidth  data  points).  Apparently,  the  energy 
requirements  for  spark  formation  were  more  time  dependent  than 
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those  required  to  increase  the  surface  scatter.  Conduction  of  thermal 

energy  from  the  irradiated  sites  during  the  duration  of  the  pulse  can 

1/2 

be  ignored  since  (D’’’  )  w  (e.  g.  ,  the  diffusivity  D  of  ZrO  is 

-  2  2  ° 

0.81  x  10  cm  /sec,  pulsewidth  is  32  nsec  and  the  spot-size  is 
72  urn). 

Now,  Bliss  and  Milam  [29]  recently  measured  the  damage 

thresholds  of  a  TiO  /SiO  multilayer  reflector  with  a  single-mode 

^  ^  2 

ruby  laser  to  be  about  2  and  50-60  J/cm  for  pulscwidths  of  20  psec 
and  20  nsec,  respectively.  It  is  reasonable  to  assume  that  the  mini¬ 
mum  threshold  of  a  single  layer  of  ZrO  for  very  short  pulses  is  also 
2 

near  2  J/cm  .  Extrapolating  the  straight  lines  in  Figure  48  to 

2 

shorter  pulses,  the  threshold  of  2  J/cm  is  reached  at  5  nsec  for 

spark  detection  and  at  1.5  nsec  for  LIS  detection.  These  pulsewidths 

are  much  larger  than  the  20  psec  used  by  Bliss  and  Milam;  so  it  is 

plausible  that  the  energy -density  threshold  could  well  be  a  constant 

value  over  the  range  of  20  psec  to  2  nsec. 

This  temporal  dependence  of  the  damage  threshold  of 

ZrO  is  similar  to  that  described  by  BHbs  [30]  for  damage  caused  by 
£ 

an  electron  avalanche  breakdown.  According  to  his  analysis,  the 
dependence  of  the  threshold  energy  density  on  pulse  duration  changes 
from  linear  with  T  for  long  pulscB  to  independent  of  T  for  short  pulBeB. 
This  result  assumes  that  the  lifetime  T  of  electrons  in  the  conduction 
band  is  long  compared  to  the  time  t‘  required  to  accelerate  an  electron 
to  the  energy  of  the  conduction  band. 
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Figure*  48.  Laser  pulsewidth  dependence  of  the  damage  threshold 
for  a  eingle  half-wave  thick  film  of  ZrC>2  on  BSC- 2 
glass.  Damage  was  detected  by  both  laser-induced 
scattering  und  spark  formation.  Deviation  from  con¬ 
stant  power  density  thresholds  is  indicated  by  dashed 
curves. 


VII.  COATING  PARAMETERS 

1.  Introduction 

Unlike  the  bulk  and  surface  damage  studies  in  trans¬ 
parent  dielectrics,  the  understanding  of  damage  in  uielcctric  thin-film 
coatings  is  exceedingly  difficult  due  to  the  variability  of  certain  film 
properties  due  to  coating  procedures.  Our  approach  to  this  multi- 
variable  problem  was  not  to  attempt  either  to  establish  any  one  coating 
procedure  or  to  control  any  particular  coating  property,  but  to  ubc 
"state-of-the-art"  coatings  obtained  from  commercial  sources.  In 
particular,  the  coating  configurations  studied  were:  (1)  single-layer 
dielectric  coatings  of  different  materials  and  thickness,  (2)  bi-layer 
antireflection  coatings  of  different  materials,  (3)  multi-layer  reflectors 
of  different  materials  and  varied  number  of  layers,  and  (4)  substrate  of 
different  materials.  In  addition,  absorption  by  thin-fi’m  systems  was 
analyzed  and  other  possible  damage  mechanisms  were  examined. 

2,  Single-layer  CoatingB 

a.  Film  Materials 

The  damage  thresholds  for  single,  quarter-wave  thick 
films  of  MgF.,,  SIC>2,  ZrC>2,  TiC>2,  and  ZnS  on  polished  BSC-2  glass 
substrates  arc  listed  in  Table  11  in  the  order  of  increasing  refractive 
index  at  6943  A.  The  thresholds  are  seen  to  monotonically  decrease 
with  increasing  film  index  (see  Fig.  49).  It  can  also  be  seen  that  the 
LIS  technique  was  the  most  sensitive  in  these  threshold  measurements 
for  6inglc-laycr  coatings.  The  correlation  with  index  is  expected, 
since  by  the  classical  theory  of  H.  Lorentz  a  low  index  is  associated 
with  a  large  energy  gap  between  the  valence  and  conduction  bands. 

This  decreases  the  likelihood  of  abrorption.  A  more  detailed  dis¬ 
cussion  on  the  role  of  absorption  as  a  damage  mechanism  can  be 
foui  d  in  Section  VII.  6. 
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T«bl*  II.  Damaga  thraahalda  for  quartar-wava  alngla  coatlnga  (nominal  pulaawldth.  11*1}  naac), 


Sampla 

Malarial* 

Rafractlra 
Indaa 
1*34)  A  I 

Spot  alau 

Peak  anargy  danalty 

by  apark  by  US 

Paak  poarar  danalty 

by  apark  by  LIS 

Patio  of 
apark 
to  US 

mm 

J/cm 

2 

OW/cm* 

01)3 

m«f2 

1.  31 

0.  033 

111-1)1 

111*1)1 

1.2-3.  ) 

1.  2-3.) 

1.0 

0123 

S102 

1.436 

0.033 

»  123 

117 

'  10.3 

10.3 

*  l.l 

O!  29 

Zr°2 

1.373 

0.  072 

61-61 

>3-44 

4.  3-4.3 

2.  4-3. 2 

1.6 

0121 

tio2 

2.  2* 

0.072 

37-71 

33-42 

4.  2-3.  7 

3.  0-3.  3 

1.7 

U20) 

ZnS 

2.  )3 

0.  072 

() 

3.3-12.  3 

4.7 

0.4-0.13 

'  3 

*AII  aubatratra  ara  BSC2  glaaa. 


Tabla  12.  Damaga  thraaholda  lor  varlooa  dim  thlcknaaaaa* 


Sampla 

Film 

Spot  alaa 

Paak  anargy  danalty 

Paak  powar  danalty 

mm 

J/cm* 

GW/cm1 

313) 

MgF2.X/4 

0.033 

111-1)1 

1.  2-3.  3 

01)3 

MgFj.k/2 

0.063 

63-102 

3. )-(.  1 

0123 

SIOj.V/4 

0.033 

117 

3.S 

0127 

SIOj.k/2 

0.033 

>  III 

<3.3 

0123 

ZrOj.  X/4 

0.  072 

36-44 

2.  4-3.  2 

01)1 

ZrOj,  X/2 

0. 072 

12-1) 

0.83-0.30 

OI2I 

TI02.V4 

0.072 

36-42 

3. 0-3.3 

0124 

TIOj,  X/2 

0. 037 

7-1 

0.33-0.63 

01)7 

TIOj,  3X/4 

0.072 

17.3-23 

1.4-2.  2 

0.  122 

OI3I 

TIOj.  3X/4 

0.072 

20-2) 

1.4-1. 3 

*Damaga  datactad  by  laaar-lnduead  tcallar. 
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INOEX  OF  REFRACTION 

Figure  49.  Correlation  of  energy  density  damage  threshold 
with  refractive  index  of  film  material. 
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b.  Film  Thickneaa 

Table  12  givea  the  damage  threaholdi  for  various  film 
thicknesses.  In  particular,  note  the  variation  of  the  thresholds  for 
the  three  TiC^  films.  The  damage  threshold  of  the  half-wave  film 
was  about  one-third  that  of  the  quarter  and  three-quarter-wave  films. 
Figure  50  calculated  relative  intensity  distribution  for  these  three  film 
thicknesses  of  TiO^.  For  a  high  index  film  material,  such  as  TiO^, 
between  two  lower  index  materials  of  air  and  glass,  the  net  electric 
field  has  an  antinode  at  the  rear  or  glass-film  interface.  This  rear 
reflected  wave  has  an  antinode  (node)  at  the  front  or  air-fiim  interface 
in  the  case  of  a  half-  (quarter)  wave  film.  The  higher  net  electric- 
field  at  the  front  surface  of  the  half-wave  film  may  be  the  origin  of 
the  lower  damage  threshold  of  the  half-wave  film.  It  is  interesting 
to  note  that  the  ratio  of  the  intensities  for  the  half  and  quarter-wave 
films  at  the  rir-film  interface  is  about  three. 

Further  understanding  of  the  thickness  dependence  may 
be  obtained  by  comparing  the  energy  densities  inside  the  films  when 
the  incident  laser  energy  was  sufficient  to  cause  threshold  damage. 

If  differences  in  the  electric  fields  in  films  of  varied  thicknesses  were 
sufficient  to  explain  the  damage  threshold  variations,  then  the  internal 
energy  densities  at  the  film  depth  most  susceptible  to  damage  should 
be  equal.  For  MgF^  and  SiO^  the  energy  densities  for  quarter  and  half¬ 
wave  thicknesses  are  most  closely  matched  at  the  maximum  values. 

For  ZrO^  and  TiO^  the  values  at  the  air-film  interfaces  are  most 
nearly  equal.  Still,  except  for  SiO^,  half-wave  energy  densities  at 
the  damage  thresholds  are  lower  than  those  of  the  quarter-wave. 
Therefore.,  even  when  electric-field  standing-wave  intensities  are  taken 
into  account,  thick  films  are  still  more  susceptible  to  damage  than 
thin  films. 


- 
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3.  Bilayer  Antireflection  Coating*  (V-coat*) 

The  damage  threshold*  for  bilayer  antireflection  coatings 

composed  of  each  possible  combination  of  the  high-index  materials* 

TiO  and  ZrO„  ,  with  the  low-index  materials,  SiO  and  MgF?,  are 
2  2  c  c 

presented  in  Table  13.  Like  single -layer  films,  LIS  was  the  most 
sensitive  method  of  detecting  damage.  The  laser  parameters  were 
essentially  identical  in  each  test,  and  the  Zr02/SiC>2  coating  was  easily 
the  most  damage  resistant.  The  damage  thresholds  of  these  coatings 
were  generally  midway  between  the  threshold  values  of  their  component 
films.  This  observation  was  also  made  by  Turner  L4]  for  coatings 
composed  of  two  quarter-wave  films.  It  is  known  that  the  fracture 
strength  of  thin  film  is  greater  than  that  of  thick  films  and  the  bulk 
material  D  3}  For  the  bilayer  coatings,  the  thicknesses  of  the  high- 
index  layers  ranged  from  only  1/5  to  1/2  of  a  quarter-wave,  and  those 
of  the  low  index  layers  were  all  about  1.  3  quarter  waves  (see  Table  6). 
This  suggested  an  explanation  to  the  fact  why  the  V-coatings  have  much 
higher  thresholds. 


4.  Multilayer  Reflection  Coatings 
a.  Film  Materials 

The  damage  thresholds  of  dielectric  reflectors  composed 
of  multiple,  quarter-wave  films  are  listed  in  Table  14.  Only  the 
photoelectrically -detected,  laser-induced  spark  thresholds  were  measured, 
because  the  initial  weak-signal  scatter  of  these  coatings  was  suffi¬ 
ciently  high  to  preclude  sensitive  LIS  detection  near  the  damage  thres¬ 
holds.  (Refer  to  Section  III.  5  for  discussion  on  weak-signal  scatter 
of  multi-layer  coatings).  The  reflectors  are  listed  by  their  film  mater¬ 
ial  combinations  beginning  with  the  most  damage-resistant  types. 
Reflectors  with  the  TiO  /SiO.  configuration  had  the  highest  thresholds, 

&  fc 

closely  followed  by  the  ZrO ^  MgF2  combination. 
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TaMa  M.  Dainap  thrraholda  for  hllayrr  a  nil  reflection  coatlnga*'** 


Simplf 

IVsr  riptlnn 

Peak  mergy  denalty 
hy  a  park  hy  LIS 

Peak  power  denalty  Hallo  of  apark  lo  LIS 

hy  aparh  hy  LIS 

J/tm^ 

J/cm^ 

GW/rm2 

OW/rm2 

014) 

7  rO^/SIO^ 

lid. 147 

47-11* 

8.6-11.7 

7. 1.. 4,  5 

1.2 

0147 

iio2/M»rJ, 

*4*5 

60.65 

5.  l-f>.  1 

4.0-4.  1 

1.4 

0145 

7.f02/M*K2 

<•4-42 

54-  57 

5.  1-6.8 

1.  5-4.0 

1.  5 

0142 

tio2/sio2 

74-48 

15-42 

f»a  5*t»a  8 

2.4-5.  1 

2.  5 

aSamplr  configuration  ia  Glll-A,  whtft 

c 

la  HSC2  (laaa  10. 5  mm  thlrkl,  11  la  high  lnilf>  film 

(<  \/4  thlcM.  L  la  low  Indei  film  1 

* 

U  4 

thlrkl,  and  A  la  air. 

^Uatr  ipo(>iltr  la  0.  Of  2 

mm,  and  nominal  pulaawtdth  la  11*14 

naac. 

Tabta  14. 

Threaholda  for 

multilayer  refleclora  (nominal  pulaewldth: 

12-15  naecl 

Sample 

Dearrlpflon 

Spot -alee 

Peak  energy  denalty 

Peak  power  denalty 

mm 

J/rm2 

GW/cm2 

0105 

TIOf/SIOf 

0.056 

121-126 

4.  1-4.  5 

0104 

0.056 

107-127 

6.  1-8.4 

*0104 

•• 

0.046 

47-118 

7. 6-8.  7 

OIOI 

•• 

0.062 

85-101 

4.  5-8.4 

0102 

0.  060 

48-110 

6.  8-7.  2 

SIOI 

0.062 

45-56 

5.0-4. 6 

SI02 

0. 062 

40 

5.7 

0107 

7-rO^/MtFj 

0.058 

81-115 

4.  7-7.8 

0108 

II 

0.062 

40 

6.  1 

0104 

F.rO^/SIO^ 

0.15 

18-20 

1.2-1. 6 

0106 

l» 

0.15 

18.  4 

«.  5 

St  01 

I* 

0.12 

7.5 

0.  52 

SSI05 

C*02/SI02 

0.  14 

17-14 

I.  4 

SSI02 

14 

0. 14 

14. 4-17 

1.  1-1.2 

SS’OI 

H 

0. 14 

11.4 

1.  I 

ssmt 

•• 

0.072 

14.4-14 

1. 0-1.4 

Ul  1 1 

7nS/ThF4 

0.084 

25-26 

2.  0-2. 1 

*  Sr  rood  laal  altar  2)  monlha. 
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As  illustrated  in  Figures  28  and  29  in  Section  IV,  the 

electric  field  intensities  peak  at  the  H-U  interfaces.  Therefore,  for 

laser  irradiation  above  threshold,  damage  involves  removal  of  lay  «*s 

in  pairs  since  the  high-index  films  are  less  damage-resistant  than  the 

low-index  films.  This  was  previously  pointed  out  by  Turner  [4]. 

For  example,  at  the  first  peak  in  the  electric  field  squared  inside  ihe 

2 

best  TiO^/SlO^  reflector,  the  energy  density  (J/cm  )  in  the  high  index 
layer  was  65%  higher  than  in  the  low-index  layer.  Taken  together 
with  the  higher  absorption  coefficient  of  the  high-index  film,  sufficient 
thermal  expansion  could  have  occurred  to  have  overcome  the  adhesive 
bond  with  resultant  film  rupture.  Reflectors  with  TiO^/ SiO^  configura¬ 
tion  had  the  highest  thresholds.  Therefore,  it  was  of  interest  to  com¬ 
pare  the  damage  resistance  of  reflectors  of  the  same  materials  but 
with  different  numbers  of  layers.  In  Table  15,  the  spark  thresholds 
of  TiC^/SlO^  reflectors  with  12  layers  are  compared  with  a  22-layer 
design.  Also  included  for  reference  are  the  thresholds  of  single¬ 
layer  and  bilayer  films  of  these  materials.  The  reflectance  of  the 
12-layer  reflectors,  98.7%,  was  not  much  less  than  99.9+%  reflectance 
of  the  22-layer  design.  The  coating  procedure,  sample  substrates, 
and  laser  parameters  were  the  same  for  the  samples  listed.  The 
thresholds  of  the  12-layer  reflectors  were  about  10  to  15%  less  than 
those  of  the  22-layer  reflectors.  Furthermore,  the  reflector  thres¬ 
holds  were  greater  than  that  of  the  single  Ti02  film  and  nearly  the 
as  single  SiO^  films. 


same 
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Table  IS.  Damage  Ihreahold*  veraua  number  ol  layer* 
of  TIOj  and  510^  film*  (nominal  pulaewtdth,  11-14  naec). 


S«  mp'i 

Film  material* 

Numb* r  of  lay  ar* 

Order  of  layer** 

Peak  energy  denalty 
by  apark  by  US 

J/cml 

J/cml 

0124 

SIOj.  J/4 

1 

GLA 

>124 

117 

0127 

SIOj.X/2 

1 

GL*A 

>111 

>111 

0121 

TIOj.  >/4 

1 

GHA 

47.71 

44-42 

0142 

TIC./SIG, 

(0.21)  11.^20)1/4 

2 

gii'l'a 

74-9* 

44-42 

0109 

T102/SI02 

12 

goim4hl2a 

94-114 

-* 

OHO 

•• 

12 

•• 

100 

-- 

C101 

•• 

22 

GmL),0IIL*A 

121-126 

•* 

0104 

M 

22 

a* 

107.127 

-  " 

0104 

•• 

22 

♦  ♦ 

97-111 

•• 

*  Ci  pottahed  DSC2  gtaaa  aubatrata,  Li  low  Indaa  film,  X/4  thick  at  6944  L  Hi  M|h  Index  film. 
1/4  thick  at  6941  L  and  Ai  air. 


5.  Damage  Threaholda  for  Coatings  on  Various  Substrates 
The  possible  relationship  of  substrate  material  with  thin 
film  damage  thresholds  was  investigated  for  quarter-wave  ZnS  films 
and  multilayer  reflectors  of  TiO^/SiO^.  The  results  are  given  in 
Table  16.  The  ZnS  films  were  depostied  on  cleaved  NaCl  crystals 
and  polished  BSC-2  glass  substrates  (samples  U201,  U202)  using  resis¬ 
tive-heating.  ZnS  films  were  also  deposited  on  BSC-2  glass  and 
polished  spinel  crystals  (U203,  U204,  U301)  by  electron-gun  evaporation. 
The  thresholds  for  the  ZnS  films  coattd  with  the  electron-gun  probably 
should  not  be  compared  with  those  coated  by  resistive  heating  since 
these  were  the  first  films  prepared  with  a  new  electron  gun  system. 

The  thresholds  for  films  on  NaCl  substrates  were  significantly  higher 
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(by  40-80%)  than  for  the  glass  substratec;  likewise,  the  threshold?  for 
films  on  spinel  substrates  had  about  60%  higher  thresholds  than  films 
on  glass.  No  significant  differences  in  the  thresholds  of  the  TiC>2/  SiC>2 
reflectors  on  polished  BSC-2  glass  and  microscope  slides  were  measured. 


Tabic  It.  Damage  thresholds  for  coatings  on  various  subalralas 


Sample 

Subelrele  description 

Spot  -  else 

Pea 

k  energy  density* 

Peak  power  density 

Single  layer  of  ZnS.k/4  thlcki 

mm 

J/cm* 

OW/cm1 

1 101 

Gfivtd  NiCl 

0.  14 

17. IT. 4 

1.3.1.  4 

U 102 

0.  164 

ll-l).  4 

1.0- 1.  1 

U20I 

Polished  f)SC2  glees 

0.  19 

9.  1.9.  4 

0.74.0.  SO 

11202 

•• 

0.  19 

9.6 

0.  SO 

1120) 

0.  14 

4.  2.4.4 

0. 10.0.40 

V  204 

• 

0.  14 

).  6-4.  1 

0.24.0.14 

•• 

•• 

0.24 

).2).S 

C. 24-0. 14 

U)0l 

Pollehed  spinel 

0.24 

4.0-6.  1 

0. 40-0.40 

Multilayer 

TIOj/SIOj  reflectors  (99«%): 

OI0) 

Polished  DSC2  glass 

0.046 

121.126 

9.  1-9.  1 

0104 

0.046 

10*. 127 

6.  1-8.9 

Ulll 

Glees  microscope  elide 

0.044 

104-121 

S.  4- 10.  S 

OII2 

0. 044 

94. 116 

7-9 

0.06) 

*  Thresholds  determined  by  laser. Induced  rcatter  (LIS)  for  single  leyere  and  by  photoelectric  detection 
of  epark  emission  for  multilayers. 


Considering  substrates  in  general,  there  are  at  lpast 
five  properties  which  can  affect  the  film  damage  threshold*  thermal 
diffusivity  D,  coefficient  of  thermal  expansion  C  ,  optical  absorptivity  8, 
film  adhesion  A  and  surface  preparation.  For  optimum  damage  resis¬ 
tance,  D  should  be  very  large,  a  should  be  about  equal  to  that  of  the 
adhering  films,  8  should  be  negligible,  A  should  be  maximum,  and  the 
density  of  defects  and  absorption  centers  due  to  surface  preparation 
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should  bo  minimal*  Cleaned  NaCl  crystals  were  expected  to  be  better 
substrates  than  BSC-2  polished  glass  since  NaCl  has  a  thermal  diffusivity 

2 

six  times  larger  than  glass  (D  =  33  x  10  cm  /see  for  NaCl  compared 
-3  2 

to  5  x  10  cm  /see  for  glass).  Both  glass  and  NaCl  have  negligible 

absorption  at  6943  A,  but  the  thermal  expansion  coefficient  of  NaCl  is 
-5  o  -1 

4  x  10  C  which  is  3  to  6  times  larger  than  that  of  several  coating 
materials.  BSC-2  glass  is  more  suitable  in  this  property  with  a  equal 
to  6.  5  x  10~^  °C**.  Adhesion  properties  of  glass  are  apparently  satis¬ 
factory,  but  the  adhesion  of  films  to  NaCl  substrates  is  subject  to  its 
hygroscopic  nature. 

The  comparable  properties  of  spinel,  a  face-centered 
cubic  crystal,  should  make  it  a  better  substrate  than  glass.  In  addition, 
it  does  not  have  the  hygroscopic  property  like  NaCl.  The  thermal 
diffusivity  is  high,  35  x  10cm^/sec  and  Q  Is  7  x  10  ^  t  *,  similar 
to  glass  and  the  films.  The  surface  of  spinel  is  not  as  easily  polished 
as  glass,  however,  as  indicated  by  the  9oA  rms  roughness  specified 
by  the  manufacturer.  A  15  to  20 A  rms  value  is  attainable  on  glass 
surfaces. 

Because  the  damage  thresholds  of  ZnS  films  on  NaCl 
were  relatively  high  despite  its  large  thermal  expansion  coefficient,  it 
is  concluded  that  either  thermal  expansion  is  not  an  important  consider¬ 
ation  for  Q-switched  laser  pulses,  or  that  the  advantages  of  a  crystal 
substrate  more  than  compensated  for  this  high  coefficient.  In  order 
to  reach  more  specific  conclusions,  further  study  must  be  conducted 
on  substrate  materials. 

6.  Evaluation  of  Absorption  as  a  Damage  Mechanism 

a.  Calculation  of  Effective  Absorption  Coefficients 

If  absorption  of  a  laser  pulse  causes  damage  by  heating 
the  film  to  its  melting  point,  what  value  of  the  linear  absorption  coef¬ 
ficient  8  1*  required?  This  question  can  be  answered  by  a  simple 
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thermal  analysis.  The  energy  absorbed  W  per  unit  volume  necessary 
to  cause  film  melting  is  given  by 

W  1  P  C  C„  -  300°K)  +  AH  ]  , 

P  mp 

where  p,  C  ,  T  ,  and  AH  are  the  film  density,  specific  heat,  melting 
p  mp 

temperature  and  latent  heat  of  fusion.  Another  expression  for  W  in 
terms  of  the  optical  field  of  the  laser  beam  is  given  by 

W(z)  =  pn 

E 

o 

where  Wq  is  assigned  the  value  of  the  axial  energy  density  at  the 
damage  threshold.  By  equating  the  two  expressions  and  rearranging 
terms,  the  absorption  coefficient  at  a  distance  z  from  the  air-film 
surface  necessary  for  meltinq  the  film  is 

--  -rfr  t  Cp  <Tmp  -  300°K> +  4»3 

o 

In  this  derivation,  conduction  of  heat  from  the  irradiated  area  was 
considered  negligible  during  the  12-nsec  laser  pi'lses,  since  (Dt)^«wq 
for  the  thin  films  tested.  Thermal  radiation  and  convection  losses 
were  also  assumed  to  be  small  during  irradiation. 

The  physical  and  thermal  properties  of  the  film  materials 
are  listed  in  Table  17.  The  density  p  of  thin  films  is  equal  to  the 
packing  fraction  f  times  the  bulk  density  p^.  The  packing  fraction  can 
be  determined  from  the  formula  [32] 

(n2-l)  (n?  +  2) 
f  =  — r - - - 

<n£-  l)(n%2) 


t 
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where  n  and  n,  arc  the  refractive  indices  of  the  film  and  the  bulk 
b 

material,  respectively. 

We  can  now  calculate  the  absorption  coefficients  that 

would  be  required  to  cause  melting  of  the  films  at  the  threshold  laser 

intensities.  Using  the  physical  parameters  given  in  Table  17,  the 

threshold  values  W  listed  in  Table  11  for  single-layer  films  and  the 

o 

relative  electric-field  intensity  distributions  given  in  Section  IV,  these 
absorption  coefficients  were  determined.  The  results  of  the  calculation 
are  given  in  Table  18.  Since  the  electric  fields  vary  with  distance  *, 
the  absorption  coefficient  required  for  melting  will  vary  with  film  depth. 

Table  17.  Physical  and  thermal  properties  of  film  materials 


Material 

Density 

Pb 

Packing 
fraction,  f 

Specific 
heat,  c 

P 

Melting 

temperature 

T 

mp 

Heat  of 
fusion, 
AH 

gm/cm^ 

J/gm-°K 

°K 

cal/gm 

MgF2 

3.2 

0.9 

1.19 

1  528 

94.7 

SiO 

2.  32 

0.93 

0.73 

2001 

35.0 

ZrO 

Ca 

5.  60 

0.89 

0.44 

2983 

168.8 

TiO, 

4.  26 

0.86 

0.71 

2093 

142.7 

Table  1 8.  Sumrru  ry  of  analysis  of  linear  film  absorption 


Material 

Calculated  3 
to  melt  firm 

Measured  0 

Published  values 
of  0 

-1 

-1 

-1 

cm 

err 

cm 

MgF2 

40-55 

<  100 

SiO 

25-30 

<  100 

<  10  [33] 

Z,o2 

200-430 

<  100 

1350  [34] 

TiO, 

160-420 

100-500 

<  40  [33] 

2 

~  1000  [35] 

ZnS 

270-1400 

10-2000  [36] 

b.  Measured  Absorption  at  6943 X 

The  calculated  values  of  p  arc  very  much  larger  than 
those  cited  in  the  optics  literature  for  bulk  crystals  in  the  visible 
spectrum  (typically  less  than  1.0  cm  *)•  To  determine  if  these  high 
values  were  inherent  film  properties,  the  transmittance  of  half-wave 
films  at  6943 A  was  measured  with  a  Cary  14  spectrophotometer. 
Within  the  resolution  of  the  instrument  (0.1  -  0.2%  in  absorption), 
no  absorption  was  apparent  for  the  SiO^,  MgF^  and  ZrO^  films. 
Minimal  absorption  was  observed  for  TiO,.  A  half-wave  film  of 
ZnS  was  not  available  at  the  time  of  measurement. 

Abeles'  formula  [33]  for  weakly  absorbing  films  was 
applied  to  these  measurements  to  compute  the  absorption  coefficient. 


n  +  n  n  i 

»  -  f  -r1-  t  -i  ] » 

n  +  n 


where  d  is  the  film  thickness  in  half-wavelength  multiples,  n  and  n# 
are  the  refractive  indices  for  the  film  and  substrate,  and  T^  and  T 
are  the  transmittances  of  the  uncoated,  nonabsorbing  substrate 


and  coated  substrate,  respectively.  Assuming  that  the  spectrophoto¬ 
meter  resolution  corresponded  to  the  maximum  possible  absorption  .or 
the  films  that  indicated  no  absorption,  3  was  calculated  to  be  less  than 
100  cm'1  for  MgF^,  SiC>2  and  ZrO^  For  the  TiC>2  film  it  was  possible 

for  3  to  be  between  100  to  500  cm-1.  The  precision  of  the  measure¬ 

ments  would  have  been  improved  if  thicker  films  had  been  available 
and  an  instrument  with  finer  resolution  had  been  used.  These  results 
can  be  compared  to  those  cf  Heitmann  [33}  who  measured  3  to  be  less 
than  10  cm'1  and  40  cm'1  at  6328  *  for  fully  oxidized  films  of  Si02 
and  Ti02.  From  previous  measurements  of  ZnS  films  by  Vlcsenko  [36], 
the  coefficient  at  6943  X  was  estimated  to  be  between  100  and  200C  cm 
depending  on  the  film  annealing  process.  He  found  that  ZnS  films, 
annealed  at  elevated  temperatures,  had  increased  absorption  for  wave¬ 
lengths  greater  than  4Gl>0  A.  This  may  explain  the  results  in  which 
the  damage  thresholds  of  ZnS  deposited  on  heated  substrates  were  much 
lower  than  those  with  no  heating. 

c.  Comments  on  Ultraviolet  Resonance  Absorption 

The  absorption  of  a  dielectric  for  visible  wavelengths  is 
due  to  the  tail  of  the  fundamental  absorption  resonance  in  the  ultra¬ 
violet  region.  The  theoretical  relationship  for  the  extinction  coefficient 
at  a  frequency  of  an  absorption  band  is  of  the  form 


.2  2.2  ,  2  2 

(vf  -  V  )  +  Vf  V 

This  relation  showr  that  nearer  the  resonance  is  to  v,  the  greater 
the  absorption  will  be.  For  ZnS  the  primary  resonance  is  at  21 50 X, 
but  there  is  a  secondary  resonance  at  3200  JL  The  resonance  of  Ti02 
occurs  in  the  near-UV  at  3200  X  ,  and  for  Zr02  it  is  at  2490*  .  For 
MgF  and  SiO^,  the  absorption  edge  is  below  2000  *  .  For  a  CeC>2 
film,  the  resonance  was  at  3000*  with  significant  absorption  extending 
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to  4600  A,  at  which  k  fell  below  0.01.  With  the  relative  damage  resis¬ 
tance  of  these  films  in  mind,  it  is  apparent  that  the  shorter  the  wave¬ 
length  for  resonance,  the  greater  the  damage  threshold.  Acutally,  this 
statement  is  equivalent  ot  the  correlation  of  threshold  with  refractive 
index  since  n  is  related  to  k. 

To  verify  that  the  thin  film  samples  used  in  the  damage 

tests  had  no  secondary  absorption  resonances  between  those  given  above 

and  the  laser  wavelength  of  6943  A,  the  absorption  was  measured  from 

3300  to  7000  A  with  the  Cary  14  spectropnotometer.  These  spectra  for 

three  film  materials  are  shown  in  Figure  51.  By  comparing  with  the 

apparent  absorbance  of  the  uncoated  glass  substrate  while  taking  into 

account  variations  in  standing-wave  reflections  and  the  index  dispersion, 

no  other  regions  of  absorption  are  discernible.  The  ZrO^  and  MgF^ 

films  are  obviously  far  from  their  resonances,  but  the  absorption  of 

TiO  increases  drastically  below  3800A  as  the  resonance  of  3200  A  is 
£ 

approached. 

To  further  demonstrate  the  presence  or  absence  of  UV 
absorption,  a  quarter-wave  thick  film  of  each  material  was  exposed 
to  the  UV  emission  (3371  A.)  of  a  pulsed  nitrogen  gas  laser.  The  laser 
emitted  a  nominal  power  of  100  kW  (~  10  nsec  pulses  at  100  pps) 
which  was  focused  on  the  samples  with  a  15  cm  lens.  No  damage 
occurred  to  the  MgF.,,  SiOz,  and  Zr02  films,  but  the  beam  vaporized 
the  TiOz  and  ZnS  films  which  have  a  large  value  of  g  at  3371  A. 

d.  Conclusions  on  Linear  Absorption 

For  some  of  the  film  materials  the  limited  precision 
of  the  g  measurements  prevents  final  conclusions;  however,  it  is 
possible  to  make  some  reasonable  remarks.  The  possibility  of  linear 
absorption  being  involved  in  laser-induced  damage  at  the  ruby  wave¬ 
length  is  greatest  for  the  higher-index  films.  This  is  almost  certain 
for  ZnS  and  likely  for  TiO^  For  MgF2  and  Si02,  however,  other 
mechanisms  such  as  nonlinear  absorption  at  the  high  intensities 
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Figure  51 
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required  for  damage  must  be  considered.  A  combination  of  both 
may  be  involved  in  the  damage  of  ZrO^  films.  To  reach  definite 
conclusions  regarding  the  role  of  linear  absorption,  measurements  of 
the  extinction  coefficient  more  accurate  than  already  obtained  ’n  this 
study  would  be  required.  It  is  not  sufficient  tc  use  values  of  0  or  k 
published  in  the  literature  because  film  properties,  including  absorp¬ 
tion,  are  affected  significantly  by  the  deposition  technique.  With  some 
techniques,  complete  oxidation  of  the  metal  ions  is  not  attained,  and 
this  deviation  from  stoichiometry  (e.  g.  Ti,  TiO,  and  Ti2°3  rat**er 
than  TiO  )  results  in  increased  absorption  at  these  "impurity"  defects. 

£m 

7.  Film  Stresses  (inherent  and  laser-induced) 

There  are  many  probable  damage  mechanisms.  Some 
of  the  results  had  characteristics  predicted  by  one  or  more  of  several 
damage  mechanisms.  For  example.  Figure  52  displays  a  damage  site 
above  spark  threshold  on  a  single  quarter-v/ave  thick  film  of  MgF^. 

The  explanation  of  this  morphology  requires  more  careful  analysis. 

It  is  possible  that  the  damage  thresholds  of  thin  films 
may  be  related  to  their  net  residual  stresses  which  are  established 
during  the  vacuum-deposition  process  and  subsequent  exposure  to  the 
atmosphere.  Ennos  [37]  and  Heitmann  [33]  have  measured  the  resi¬ 
dual  stresses  in  a  number  of  coatings  including  TiO^,  SiO^,  MgF^, 
and  ZnS.  The  stresses  in  TiC>2  and  MgF£  are  tensile  with  magnitudes 
of  3  x  10^  and  1  x  10^  dynes/cm2,  respectively.  Those  of  SiO  and 
ZnS  are  compressive  with  magnitudes  of  about  2  x  10  dyneu/cm  . 

The  stresses  in  most  other  materials,  including  Zr02  are  tehsile. 
Stress  compensation  may  be  effected  by  proper  pairing  of  ".Ims  with 
compressive  and  tensile  stresses.  Since  very  thick  films  have  been 
observed  to  cloud  or  crack  without  applying  other  external  forces  [38  , 
it  is  apparent  that,  without  compensation,  the  magnitudes  of  these 
stresses  can  be  of  the  order  of  the  fracture  strengths.  As  a  result, 
the  magnitude  of  additional  forces  required  to  cause  film  rupture  may 


Figure  52.  Damage  site  on  a  single  quarter-wave  film  of  Mg 
with  (a)  overall  view  of  site  having  dimension  of 
150  microns  diameter  and  (b)  close-up  view  of  the 
far  right  region  of  the  damage  site. 
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be  substantially  modified. 

Sinci  tht  radial  distribution  of  a  laser  beam  is  necessar¬ 
ily  nonuniform,  e.  g  Gaussian,  thermostrictlve  stresses  result  from 
film  absorption.  The  temperature  rise  in  the  thin  film  due  to  absorp¬ 
tion  of  a  Gaussian  beam  is  given  by 

,  2,  2 

AT(r)  =  AT  e'2r  /w  , 
o 

and  the  resultant  stress  components  for  steady-state  absorption  are 
given  by  [39] 

«r  ■  •°Y4To(ir]2  C>  -e'2r2/"Y 

and 

,2.2 

°e  ’  -°Y4To'  /w  -  ;r 

where  o  is  the  coefficient  of  linear  expansion,  and  Y  is  Young's 

modulus.  The  radial  stresses  are  compressive,  whereas  the  azimuthal 

stresses  are  tensile.  Although  o  is  always  larger  than  o  »  the  com- 

0 

pressive  strength  of  materials  >s  much  greater  than  the  tensile  strength 
(by  a  factor  of  10  or  more  for  glass,  for  example)  [38], 

Using  a  dynamic  analysis  of  thermostriction,  Sharma 
and  Rieckhoff  [40]  determined  that  mechanical  damage  of  the  bulk  of 
silicate  glasses  whould  result  from  tensile  stresses  if  the  effective 
absorption  coefficient  was  at  least  50  cm”*.  (Surface  damage  would 
require  a  smaller  coefficient).  Recalling  the  large  absorption 
coefficients  listed  in  Table  18  postulated  for  film  melting,  it  is 
reasonable  that  thermostrictive  forces  could  also  exceed  the  fracture 
strength  of  thin  films  with  those  values. 

Electric  surface  stress,  electrostriction,  and  radiation 
pressure  are  other  laser-induced  stresses  which  have  been  considered 
as  possible  causes  of  surface  damage  [41],  Using  the  theory  of 
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Stratton  [42],  Kerr  T41]  derived  an  expression  for  the  electric  surface 
stress  as 


s  =  <nc  +  no  P12  '  Uno)I/2C  ' 

where  S  is  the  pressure  directed  outwardly  normal  to  the  surface  PJ2 
is  of  the  order  of  0.25,  I  is  the  laser  intensity  (watts/cm  ).  and  c  is 
the  velocity  of  light.  For  a  laser  pulse  with  peak  intensity  of 
10  GW/cm2  (the  maximum  value  measured  for  the  most  damage  - 
resistant  coating),  the  electric  stress  is  equal  to  5  atmospheres. 

This  is  quite  small  when  compared  to  the  value  required  to  damage 
glass  (-100  atm)  [40].  Radiation  pressure  and  electrostriction  within 

a  film  would  be  even  lower. 


VIII.  SUMMARY  AND  CONCLUSIONS 


Using  a  single  transverse  and  longitudinal  TEMQ0, 

Q  switched  ruby  laser,  damage  phenomena  have  been  studied  in 
dielectric  thin-film  sy  '♦errs.  This  study  included  single-,  bi-, 
and  multi-layered  vacuum-deposited  coatings  of  the  materials  MgF  , 

Ct 

SiO^,  ZrO ,  TiO^,  ZnS  and  others  on  substrates  of  glass,  fused 
silica,  rock  salt,  ana  roinel.  Damage  thresholds  were  measured 
using  spark  detection  system,  LIS  technique  and  microscope  as  a 
function  of  many  parameters  of  the  laser  beam  and  the  coatings. 
Major  experimental  results  are  summarised  in  the  following: 

1.  The  damage  thresholds  of  thin-film  systems,  as 
detected  by  increased  scattering  of  a  He-Ne  gas  laser  (LIS)  occurred 
before  or  at  the  spark  formation  detectable  by  sensitive  photoelec¬ 
tronics.  The  presence  of  a  spa  k  gives  a  very  characteristic  and 
unmistakable  damage  morphology  on  the  coatings.  The  LIS  tech¬ 
nique  was  most  sensitive  for  detecting  damage  on  single-  and  bi- 
layered  coatings.  Using  the  indication  of  the  weak-signal  scatter 
level  of  an  attenuated  ruby  beam,  the  damage  thresholds  were 
generally  higher  for  multilayer  reflectors  with  the  least  scattering. 
Time-resolved  measurements  of  scatter  (a  He-Ne  probe  beam)  have 
shown  that  the  thin  film  damage  occurred  within  the  ruby  pulsewidth 
of  20  nsec. 


2.  The  damage  thresholds  of  thin  films  were  strongly 
dependent  on  the  standing-wave  patterns  of  the  internal  electric  fields. 
The  entrance-face  thresholds  were  equal  to  or  greater  than  the  exit- 
face  thresholds  of  thin  films.  Calculation  of  the  electric-field  dis¬ 
tribution  in  thin-film  systems  was  essential  to  the  analysis  of 
damage  processes. 
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3.  In  correlating  the  damage  thresholds  to  laser 
parameters,  one  of  the  important  results  cf  this  study  is  that  the 
damage  threshold  of  the  thin  film  increases  as  the  spot-size  of  the 
laser  beam  decreases.  A  simple  model  has  been  developed  in 
correlating  the  nature  and  distribution  of  coating  defects  to  this  spot- 
size  dependence.  It  has  been  demonstrated  that  this  is  a  good  model: 
the  probability  of  the  laser  beam  striking  a  defect  site  is  greater  for 
larger  spot-sizes  while  damage  in  materials  can  be  distinguished  as 
defect  damage  and  intrinsic  damage. 

4.  The  damage  thresholds  of  coatings  irradiated  by 

a  multimode  laser  were  significantly  less  than  by  a  single-mode  laser. 
For  10  to  30  nsec  laser  pulses,  energy  density  rather  than  power 
density  was  the  more  proper  measure  of  the  damage  threshold.  The 
energy-density  threshold  for  a  ZrC>2  film  increased  with  increasing 
pulse  width,  but  less  than  that  giving  a  constant  power-density  threshold 

5.  For  single-layer  dielectric  coatings,  the  damage 
thresholds  were  inversely  proportional  to  refractive  indices  of  the 
materials.  MgF2  had  the  highest  damage  resistance  among  all 
tested  samples.  The  damage  thresholds  of  half-wave  thick  films, 
except  Si02,  were  less  than  those  of  quarter-wave  films,  in  good 
agreement  with  electric  field  predictions  at  the  air-film  interface. 

6.  Thresholds  of  multi-layer  coatings  with  high  damage 
resistance  were  much  larger  than  the  thresholds  of  the  low-index 
component  films  and  less  than  or  equal  to  the  thresholds  of  the  high- 
index  components. 


7.  Damage  thresholds  of  single-layer  ZnS  films  on 
crystalline  substrates  (cleaved  NaCl  and  polished  spinel)  were  greater 
than  on  polished  glass  substrates. 
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The  analysis  indicated  that  significant  linear  absorp¬ 
tion  could  raise  the  film  temperature  to  the  melting  point  or  high 
enough  for  thermostrictive  forces  to  exceed  the  tensile  fracture 
strength  of  the  thin  films.  For  high-index  films,  such  as  ZnS. 
linear  absorption  appeared  to  be  the  probable  mechanism.  However, 
it  was  not  evident  that  the  thresholds  of  low-lnde<  films,  such  as 
MgF2  or  Si02,  were  established  by  linear  absorption.  The  spot-size 
effect  suggested  that  local  coating  defects  are  playing  an  important 
role  in  the  damage  processes.  Using  the  defect  model,  further 
studies  on  the  parameter  of  enhancement  factor  T| ,  and  hence  either 
the  intrinsic  damage  threshold  Ij  or  the  defect  damage  threshold  , 
of  the  materials  will  uncover  one  or  more  of  the  damage  mechanisms 
in  the  dielectric  thin-film  systems.  For  other  laser  wavelengths  of 
interest,  such  as  1.06nm  and  10.6^m,  the  role  of  local  defects  via 
the  spot- size  dependence  should  be  explored. 


APPENDIX  A 


CALCULATION  OF  ELECTRIC  FIELDS  IN 
THIN  FILM  COATINGS 

The  configuration  A  the  thin  film  layers,  and  the 
notation  to  be  used  in  calculation  of  the  electric  fields  are  presented 
in  Figure  53.  As  shown,  the  fields  and  interfaces  are  numbered 
starting  from  the  surface  nearest  the  incident  wave.  For  the  con¬ 
venience  of  the  readers,  two  computer  programs  (in  BASIC  language) 
are  included  at  the  end  of  the  Appendix.  The  square  of  the  absolute 
values  of  the  electric  fields  for  different  cases  are  summarized  in 
the  following: 

1.  Single  Layer  Coatings 
a.  Normal  Incidence 

IE  /  E+  | 2  =  (A2+  B2)  /  |Dj  2  ,  (A.l) 

1  o  o' 

jEj/ E^|2  =  t2  [1+ r2  +  2r2  cos2  (bj-kjz)  ]  /  |  D  |2  ,  (A.  2) 

»nd  |E2/E*|2  =  «2.2  /  |D|2  .  (A.3) 

where  A  =  ( l+r2>  [cos(6j  -  kQz)  +  ^  cos  (6j  +  kQz)  ]  , 

B  =  (l-r2)  [  sin  (6j  -  kQz)  +  rJ  sin  (6j+  kQz)  ]  , 

|  D|2  =  1  +  r2  r2  +  2rj  r2  cos  26j  . 
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b.  Non-normal  Incidence 

For  example,  p-polarization  is  given  here 
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IE  /E+|2  =(A2+B2)  /  |  D| 2  .  (A-4) 

1  o  o' 

lEl/ Eo^2  S  C1  +  r2  +  2r2  COS  <2kl*  ’  26J  +n)  ^  /  nl  lD'2  (A,5) 

and  |E2/E^,2  =  t^t2  /  n2,D,2  .  (A.  6, 


where  A  =  <1  +  r2>  f  cos<  ^  -  kQ  z-tt  )  +  r }  c03(6j  +kQZ  +rt  )  ]  , 

B  =  (1-  r  }  f  sinlftj  -  kQz  - tt )  +  Tj  sini&j  +  kQz+n)  ]  , 
2  2  2 

and  |D|  =  1  +  ^  r2  +  Zr}r2  cos  2  6j  . 


2.  Bilayer  Antireflection  Coatings 
In  the  first  layer,  we  obtain 

|E.  /  E+  1 2  =  (G2+H2)  /(E2  +  F2)  ,  (A.  9) 

i  1  o' 

where 

G  =  tj  [coskjZ  +  r2r3cos  (262+kjZ)  +  r2cos(26,  -  k^z) 

^  ^  jC08(26j  4*  2^2^]  z)  3  9 
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H  =  tj  T sin  kjZ  +  r?r?  sin(262  +  kjZ)  +  r2  sin(26j  -  kjZ) 

+  r ^  sin(26j  +  26  2  -  k^z)  ]  . 

E  =  1  +  rr  cos  2*>j  +  rjr3  cos  2(6j  +  62)  +  r2r3  sin  2&2  . 

F  =  r j r2  si  .  26.(  +  rJr;}  sin  2< 6j  +  62)  +  sin  2&2  , 

The  respective  fielda  in  the  c?.-ond  layer  and  in  the  substrate  are 

IE  /E+|2  =  t2  t2  (1+r2  +  2r  cos  ?-(6--kz)  /(E2+F2)  (a.  10) 

1  2  O  1  12  3  3  22 

and 

|E  /  E+ 12  =  t2  t2  t?  /  (  EZ+F2)  .  (A.  11) 

’  3  o  !  12  3 

3.  Multilayer  Reflectors 

Electric  field  distributions  for  multilayer  reflectors 

were  computed  by  the  matrix  method  for  normal  incidence.  Such 

reflectors  are  usually  composed  of  alternating  layers  of  two  materials 

with  different  refractive  indices.  For  quarter- wavelength  thick  films, 

nd  =X/4,*o6  =kd  =  l/2n,  except  6  =  0  as  before.  The 

mmo  mmm  o 

matrix  relation  for  the  electric  field  wave  amplitudes  at  the  boundary 
of  the  (m-l)th  and  the  mth  layers  is 


(A.  12) 
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where  0  <  m  <  M+l  and  M  is  the  total  number  of  layers.  For  the 


incident  medium,  m=0 ,  bo  &q  =  0  and 


K  3  C 


(A.  13) 


For  the  substrate,  m  =  M+l,  in  which  there  is  no  reflected  wave, 
i.  e. .  E^+1  =  0,  ..  that 


EM  *  (l/tM+lIEM+l 


(A.  14) 


(irM+l  /  *M+1  *  EM+1 


(A.  15) 


To  evaluate  E*  including  E+  by  computer,  it  was  necessary  to  set 
M  o 

E+  equal  to  a  constant  number.  When  the  resultant  numerical 
M+ 1  +  +  + 

value  of  the  incident  field  Eq  was  determined,  Em  and  em+i  were 

then  properly  normalized  by  division  by  E*  .  As  in  the  preceding 

cases,  the  electric  fields  within  each  layer  are  giver,  by 


_+  -ik  z  ,  +ik  z 

E  =  E  e  m+Ee  m 
mm  m 


which,  for  the  present  case,  yields 


I-ik  z  +ik  z 

,  m  m 

Em+1  «'  "rm+l  *  ' 

-  ik  z  +ik  z  *| 

.  m  m 

+  Em+1  rm+1  '  -e  ’  I 


(A.  16) 
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If  a  half-wave  thick  film  is  used  for  the  first  layer,  which  is  a 
common  practice,  it  may  be  treated  as  two  adjacent  quarter-wave 
thick  films. 


Program  1  HULFIL:  Calculation  of  the  Electric  Field 
Intensity  Distribution  for  Light  Incident  Upon  a 
Multilayer  Coating  of  N  Quarter-Wave  Thick  Films 

LOAD  (HULFIL) 


7LIST 

10,  DECLARE  N  CONTROLLED  ■>  A  CONTROLLED  .  D  CONTROLLED  1 

?0.  IF  ALLOCACN)  THEN  FREE  N 1 

30.  IF  ALLOCA(A)  THEN  FREE  A) 

AO.  IF  ALLOCA(D)  THEN  FREE  Bl 

50*  CET  L1STCM)! 

60.  ALLOCATE  NCM). A(M).DCM)! 

70.  GET  LI ST(N> I 

80.  LET  RCI)eCNCI-l)-N(I))/(NCI-l)+N(I>>! 

90*  BCM>»01 

100*  A(M>»1 J 

HO.  LET  CCI)-ACI)/AC1)1 

180.  LET  DC  I  )*>B(  I  )  /A(  1>I 

130*  L00P1 t  DO  1*1  TO  M-2! 

1  AO*  ACM- 1 )*CACM+1-I )+RCK+l-I)+BCM*l-I>)/C l+RCM+l-I  )>1 

150.  BCM-I  )«*-(  ACM+l-I  >*RCM+1-I  XBCM+l-I  >>/C  1  +  RCK+.1  - 1  >  >  * 

160*  END  I 

170.  ACl)«CACft>+nC2)*BC2>>/Cl*nC2>>I 

180*  BC  1)-CAC2)*RC2>  +  B(2)>/(1  +  R(2>>! 

190*  LET  E$1CI)*CCCI)+DCI))**21 

800*  LET  FS2C I)«CCI)»*2+DCl>**2*2*CCI>*DCI >*COSDC  30)1 

210.  LET  ES3C  I  )«*CC  I  )**2+DC  I  )**2+2*CC  I  )*DC  I  >»COSDC  60>* 

280.  LET  ESAC I )*C( I )**2+C( I >**2! 

230.  LET  ES5C I )*C( I )**2+D(I >**2+2*C( I >+DC I >*C05DC 120)1 

2A0.  LET  ES6CI)-CCI>**2*DC!>**2*2»CCI>»DCI>*C0SDC 150)1 

250.  LET  ES7CI)-CCCI)-DC1))**2) 

260.  PUT  LI STC  * ' )1 

270.  PUT  LI  STC*  E.  FIELD  SQUARE  IN  I  -  T>I  MEDIUM')* 

2B0.  PUTLISTC'*)1 

290.  PUT  EDITC  1.E57C  1  >.  ES6C  1  >#ES5C  1  >»  ESAC  1  >.  ES3C  1  >.  ES2C  1  >.  ES 

C  l  )  >CF(2).F(F.3).F(8.  3). BCR. 3).  FC  8. 3)#  FC8.  3>.F( 6# 3>. FC 8. 3) > I 
300.  L00P2*  DO  1*2  TO  K-ll 

310.  PUT  EDITC  I.ES1CI).  ES2C  I  ).ES3(  I  )r>ESA(  I  >  j  ES5C  Dr  ES6C  I  )#  ES 

C I ) MFC  2). FC8.3).FC8. 3).FC8. 3) . FC S*^) . F( 8. 3># FC8. 3)> FC 8* 3>>* 

320*  END  1 

330.  PUT  EDITCM.ES1CM)  )  C  FC  2)#  FC8*  3)  )  1 

3A0*  LET  NES1CI  >«N( I) *ES1CI >1 

350.  ET  NES2C I > *N( I >*ES2( I )1 

360.  LET  NES3C  I)*N(I)*ES3(I>1 

370.  LET  NESAC  I  )*NC  I )  +  ESAC  1)1 

380.  LET  NES5C 1 )*N( I >+ES5C 1)1 

390.  LET  NES6CI)»NCI)*ES6CI)1 

AGO*  LET  NES7C I >“NC I )*ES7C 3 )1 


CONTROLLED 


CONTROLLED 


L00P1 * 


330* 

3A0* 

350. 

360. 

370. 

380* 

390* 

AGO. 


AIO.  PUT  LI  STC  f<  *  >1 

A20.  PUT  LI  STC  *  N*E.  FIELD  SQUARE  IN  I-TH  MEDIUM')! 

A30.  PUT  LIST*")! 

AAO*  PUT  EDI  TC  1  .NES7C  1  >  » NES6C  1 ) . NES5C 1 ).  NESAC  1  ).  NES3C  I  ).NES2 

1  >.NES1  (l))(F(2>.F(8.3).)Cjn.3>.F(8.3>.F(8*3>.FC8.3)/F(B.  3>#F(8*3)>! 

A 50*  LOOPS!  DO  1*2  TO  M-ll 

A60  •  PUT  EDI7  C  I  .NES1  CD#  NES2C  I  )  .NE53C  I  )  .NESAC  I  >.  NES5C  I  ).  NES6 

I  ).NES7(  I))(F(2).F(R.  3)»F(8.3).F(8.3I.F(8»3>.F(8.3).F(8.3).F(8.3>)1 
A70*  END  1 

ABO*  PUT  EDITCM.NESICK))CFC2).FC8.3>>! 
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110. 


Program  2  VCOAT:  Calculation  of  the  Electric 

Field  Intensity  Distribution  fer  Light  Incident 
Upon  a  Bilayer  Antireflection  Coating 


LOAD  C VCOAT) 

7LIST  10  THRU  246 

10.  GET  LISTCNO.Nl  #N2#  N3»  01 . 02»  D1  #  D2)  ) 

20.  Rt-CNO-Nl  )/CNjO+N1  )) 

30.  R2=CN|-N2)/<N1+N2)I 

AO.  n3«C:i2-N3)/<N2+N3)l 

50.  Tl-l+Rl) 

60*  T2°>  *R2J 

70*  T3^'-'A31 

00.  E«l4Ri*R2*C0SD(?*Cl)+Rl  *n3+C05D( 2+01+2+C2)+R2*R3+C0SDC2 


F"RI*R2*SIND(2*Q1  )  +  RH-R3+SIND<2*Ql+2+02)+R2+R3+SINDC2*Q 
2 

D-E++2+F++2I 

LET  EFS0(Z)«CACZ)**2+E<Z)**2)/D) 

LET  A<Z)»( 1+RI  )*C0SD(K0+Z)+R1+F.2*C0SDC2*QI-K0*Z)  +  R1*R3* 
0SD(2+Ql+2+Q2“K0*Z)+P2*R3*C0SD<  2+G2-C0+Z  )  +  R2+COSD(2*Gl+KO*Z)+R3*COSB(  2*0 
♦2*Q2+K0*Z)+RI  *R2*R3*C0SIV2*02+K0*Z  )  I 
130.  LET  BiZ)"<Rl-l  )  +  SIND(K0+Z)+Rl*RP*S!NDC2*0i-K0*Z)  +  Rl*R3* 

IND(2*OI+2*G2-KO+Z)+R2«3i3*SIND<2*Q2-KO*Z>  +  F2*SlNDC  2+0 1  ♦  K0*Z)«‘R3*SINDC  2*Q 
♦2+Q2+K0*Z>+Rl+R2*R3*SlNL<2+Q2+K0*Zp; 

1A0.  LET  EFS1  CZ)«CG<Z)**2+H<Z)+*2)/PJ 

150.  LET  G(Z)cTl*  CCOSDCKI  +  Z)  +  R2+K3*C0SD( 2*Qi-.+Kl*Z)  +  n2+C0SD(2 

0  ]  »K1  +  Z  )+R3+C0SDf  2+01  +  2+02'*Kl+Z)  )  *  ♦ 

150.  LET  HCZ)*T1  *(  SlND(Kl*Z)+R2*R3*5IND<2*G2+Kl*Z)  +  n2* SINDC2 

Ql-Kl*Z)+R3*SIND<2*CH-2+G2-Kl+Z4)J 

170.  LET  EFS2CZ)«<TI  +  T2)  ++2*<l  +  R3*-+2+P+n3*C0SDC2*02-2*K2*Z)> 

D  /  / 

|80.  EFS3-CTl  +  T2*T3)  ++2/D» 


02)  > 
90. 
)  I 

100* 
1 10* 


190. 

L-. 69431 

200. 

K0«360*N0/LI 

210. 

KI-360+N1/L1 

220. 

K2*360*N2/LI 

230. 

K3-360+N3/L) 

235. 

PUT  LISTC  *  * )I 

238. 

PUT  IMAGE(E.F#D)CIM1>1 

240.  IMlt 

IMAGE) 

242. 

PUT  LI STC '  *)) 

244. 

PUT  L1STC  *  Z 

FIELD  SQ*  ' ) 1 

E 

246. 

PUT  LISTC * 

D  * 

N+Z  E.  FIELD  SQ.  N* 

MEDIUM  NO »)) 


/ 
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111. 


VCOAI — Continued 


AA-ACZJ1 
BB-BCZ)! 

ESO-EFSOCZJJ 

NOF.SO-NO*EFSOCZ>I 

N0Z-N0*Z*4/LI  „  , 

PUT  EDIK  Z»  NOZ.  ESO#  HOESO  )  C  FC  6i4)»X(6>#F(6#4)#XUl)>F(  6» 


LIST  250  THRU  580 

250.  LOOP!*  DO  Z-0  TO  *347 1  S' BY  .0347 IS ^ 

260*  AA-A'T'»  no 

270* 

280. 

290. 

300. 

310. 

>#XC12>#FC6#4>>I 
320. 

325. 

328. 

330*  L00P2* 

340* 

350. 

360. 

370* 

380* 

390. 

) #XC 12>.FC6.4>>! 

400. 

405. 

408. 

410*  L00P3I 
420. 

430. 

440. 

“450. 

)  #XC 12>#FC6.4>>! 

460. 


MEDIUM  Nl')J 


END  I 

PUT  LI STC *  ')» 

PUT  LI STC • 

DO  Z-0  TO  Di  BY  Dl/101 
GG-GCZ)! 

HH-HCZ) ! 

ES1-KFSICZ) I 

N1ES1»N1*EFS1CZ)J 

NIZ-N1+Z*4/LJ 

PUT  EDI  TCZ.N1Z#  ESl.Nl  ESI  )(F(6#4)#X<6)#F<6#4)#XCll>#F(  6. 

4 

END  I 

PUT  LI STC *  *>J 

PUT  LI  STC*  MEDIUM  N2*)I 

CO  Z-0  TO  D2  BY  D2/5I 
ES2-EFS2CZ)  I 
N2ES2-N2*EFS2CZ)I 

N2Z-N2*Z*4/LI  %  * 

PUT  EDI  TCZ.N2Z#  ES2.N2ES2)  C  FC6.  4>#XC6)#FC6#4)#XCI  l)#FC  6. 


470. 

480* 

490. 

500. 

S10* 

520. 

530. 

540. 

550. 

560. 

570. 

580. 


MEDIUM  N3’>! 


IM2t 


END  J 

PUT  LI STC*  *>* 

PUT  LI STC  * 

N3KFS3-N3*EFS3J  ^  ^  .... 

PUT  EDI  TC  EFS3*  N3EFS3) CXC29J.FC6#  4>»XC  12># FC  6#  4>  >1 
PUT  LI STC*  *)l 
T-N3EFS3I 

PUT  LI STC  *  TRANSMI TTANCE.T*  >1 

PUT  IMAGEC T> C IM2) 1 

IMAGE! 

X  - 

PUT  LI STC*  *>! 

PUT  LI STC  *  *>! 

END  I 
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